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^ Abstract 

Coulomb dissociation is an especially simple and important reaction mechanism. Since the 
■^j- \ perturbation due to the electric field of the (target) nucleus is exactly known, firm conclusions can 

be drawn from such measurements. Electromagnetic matrixelements and astrophysical S'-factors 
for radiative capture processes can be extracted from experiments. We describe the basic elements 
of the theory of nonrelativistic and relativistic electromagnetic excitation with heavy ions. This is 
contrasted to electromagnetic excitation with leptons (electrons), with their small electric charge 
and the absence of strong interactions. We discuss various approaches to the study of higher order 
electromagnetic effects and how these effects depend on the basic parameters of the experiment. 
The dissociation of neutron halo nuclei is studied in a zero range model using analytical methods. 
We also review ways how to treat nuclear interactions, show their characteristics and how to avoid 
them (as far as possible). We review the experimental results from a theoretical point of view. Of 
special interest for nuclear structure physics is the appearence of low lying electric dipole strength 
in neutron rich nuclei. Applications of Coulomb dissociation to some selected radiative capture 
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reactions relevant for nuclear astrophysics are discussed. The Coulomb dissociation of S B is relevant 
for the solar neutrino problem. The potential of the method especially for future investigations of 
(medium) heavy exotic nuclei for nuclear structure and astrophysics is explored. We conclude that 
the Coulomb dissociation mechanism is theoretically well understood, the potential difficulties are 
identified and can be taken care of. Many interesting experiments have been done in this field and 
many more are expected in the future. 
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1 Introduction 



Nuclear reactions at high energy are very complicated in general, due to the strong interactions between 
the colliding hadrons. However, in the special case of very peripheral collisions the nuclei do not touch 
each other and the short range strong forces between the collision partners are avoided. The reaction 
mechanism becomes simple and unambiguous conclusions can be drawn from such experiments. In very 
peripheral collisions the nuclei interact with each other through the time-dependent electromagnetic 
field caused by the moving nuclei. Especially for heavy nuclei these fields are very strong and various 
interesting effects can occur. It is the purpose of this review to describe them. We discuss how these 
collisions are treated theoretically and how they are applied to nuclear structure and nuclear astrophysics 
problems. Coulomb excitation has been a very powerful tool in the past to study electromagnetic 
matrix-elements in nuclei. Classical review papers exist, see, e.g., P3I2]- For collision energies below 
the Coulomb barrier the condition of no nuclear interactions of the nuclei with each other is very well 
fulfilled and valuable nuclear structure information has been obtained. An important example is the 
investigation of nuclear rotational and vibrational collective motion by means of Coulomb excitation 
PP. Such studies are also of special interest for the investigation of exotic nuclei at the radioactive beam 
facilities which have become available in the past decades all over the world. Nuclei (far) away from 
the valley of stability (neutron- or proton-rich) can now be investigated experimentally. They have 
widened considerably the landscape of nuclear physics and also made possible novel studies in nuclear 
astrophysics. New phenomema like halo nuclei have been discovered. Electromagnetic excitation and 
dissociation is again a very powerful tool in this case. They are expected to play an even more important 
role in the future at the new facilities being proposed presently around the world (GSI, Germany [3], 
RIA, USA [HE], RIKEN, Japan 0). 

The condition that the hadrons do not touch each other can be fulfilled either by using bombarding 
energies well below the Coulomb barrier, or by going to very forward scattering angles in intermediate, 
(v < c, 7 > 1) and relativistic ( v ~ c, 7 ^> 1) collisions. These forward angles (classically) correspond 
to impact parameters larger than the sum of the nuclear radii for high beam energies. The two cases 
are shown schematically in Figs. ^ (a) and (b). 

It is very important to note that with increasing beam energy higher lying states can be excited with 
the Coulomb excitation mechanism. This can lead to Coulomb dissociation, in addition to Coulomb 
excitation of particle-bound states. This was reviewed some time ago in [Jj. It has become more 
and more clear, that such investigations are also well suited for secondary (radioactive) beams. An 
(unstable) fast projectile nucleus can interact with a high Z target nucleus. In this way the interaction 
of an unstable particle with a (quasireal or equivalent) photon can be studied. A similar method is used 
in particle physics, where it is known as the Primakoff effect [HJ E] ■ 

Since the electric field of a nucleus with high charge number Z is much stronger than, e.g., the one 
of an electron, the nucleus can be a very suitable electromagnetic probe for certain cases. One can 
study, e.g., higher order phenomena, which are inaccessible with conventional electromagnetic probes 
like the electron. The excitation of the double phonon giant dipole resonance observed at GSI [TUJ HI] 
is an example. 

We now give a guided tour through the content of this review: we start with a general discussion 
of the theory of Coulomb excitation and dissociation. Due to the time-dependent electromagnetic 
field the projectile is excited to a bound or continuum state, which can subsequently decay. Since the 
electromagnetic interaction is very well known, rather precise unambiguous information can be obtained 
from such experiments. If first order electromagnetic excitation is the dominant effect, experiments can 
directly be interpreted in terms of electromagnetic matrix elements which also enter, e.g., in radiative 
capture cross-sections. It is therefore very important to know the conditions of validity of the first order 
theory and how to assess quantitatively the higher order effects. 

We especially mention two aspects of multiple electromagnetic excitation: it is a way to excite new 
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(a) 




(b) 





Figure 1: (a) A projectile with charge Z\e is scattered on a target nucleus. The projectile velocity 
is nonrelativistic (v <C c). The minimum distance of closest aproach of the Rutherford trajectory is 
r min > R\ + R2, the sum of the nuclear radii, (b) A projectile with velocity v ^ c passes by a target 
nucleus with an impact parameter b > R\ + R2. The Rutherford orbit is very close to a straight line. 



nuclear states, like the double phonon giant dipole resonance [T2] ; but it can also be a correction to the 
one-photon excitation ^3JEJE3- The theoretical methods to treat these problems are discussed in 
Sec. El They range from higher order perturbation theory to Glauber calculations, CDCC (Continuum 
Discretized Coupled Channel) and solving numerically the time-dependent Schrodinger equation. 

Breakup reactions are at least a three-body problem, which is very complicated in general, and 
the long range character of the Coulomb interaction adds to the complication. We discuss here the 
Coulomb dissociation of neutron halo nuclei as an illustrative example. We can study analytically a 
rather simple model of the dissociation of a two-particle system (deuteron, neutron halo nucleus) bound 
by a zero-range force in the nuclear Coulomb field. This sheds light on the conditions of validity of the 
usual treatment of the semiclassical Coulomb excitation theory and of the importance of higher order 
effects. This is done in Sec. HI 

Coulomb excitation is a very useful tool to determine nuclear electromagnetic matrix-elements. 
This is of interest for nuclear structure and nuclear astrophysics [THl H3 [TBI H2] • In the following we 
illustrate the theoretical concepts and discuss examples of experimental results and their theoretical 
interpretation. The emphasis is more on theory. As theorists we cannot give a detailed treatment of 
nor do full justice to the ingenious experimental developments which make possible the precision studies 
at these extreme forward angles. We refer to the previous reviews, see e.g. |17| ITS] and to the new 
proposals for future facilities like [3] or jl]. 

In Sec. |H1 we review the results that have been obtained on nuclear structure from intermediate 



energy electromagnetic excitation in the past years. We briefly mention the Primakoff effect in particle 
physics, where the nuclear Coulomb field is used to study the interaction of fast (unstable) particles with 
the quasireal photons due to the Coulomb field of a nucleus. We also mention the very large effects of 
electromagnetic excitation in relativistic heavy ion collisions in fixed target experiments. These effects 
become quite spectacular at the heavy ion colliders RHIC and LHC with the very high values of the 
Lorentz parameter. Recently experimental results from RHIC have become available. 

In the next section we discuss applications to radiative capture reactions (the time reversed process 
of photodissociation) of astrophysical interest. This subject has been reviewed before ^H]- This field 
is rapidly expanding; this is essentially due to the new experimental possibilities world-wide. Future 
applications are given in Sec. |H1 We discuss new possibilities, like the investigation of r- and rp-process 
nuclei, which will be produced with high intensity at the future RIA facilities, and the experimental 
study of two-particle capture processes by means of Coulomb dissociation. We close with conclusions 
and an outlook in Sec. El This review grew out of a series of lectures given by one of us in the frame of 
the "European Graduate School Basel-Tubingen" , for the transparencies see [TH] . 

2 Theory of Electromagnetic Dissociation 

First we briefly recall the basic points of the theory of nonrelativistic (NR) Coulomb excitation. Ex- 
cellent reviews exist, see e.g. [T]. It is usually a very good approximation to treat the relative motion 
between projectile and target classically (semiclassical approximation) . It is also well understood how 
one can obtain the semiclassical limit from the quantal theory j2J. We want to concentrate here on 
intermediate and relativistic energies. In these cases one can replace the Rutherford trajectory by a 
straight line to a good approximation. In contrast to low energy Coulomb excitation retardation is now 
an important effect and the long wave-length approximation is not necessarily a good approximation. 
Deviations from the straight line approximation have to be and can be assessed quantitatively. In 
collisions above the Coulomb barrier strong interactions between the nuclei cannot be avoided. How- 
ever at very forward angles, corresponding to large impact parameters, such nuclear effects can become 
negligible. Even then there can be some kind of nuclear effects at a certain level of accuracy: there is 
diffraction due to the wave nature of the projectile. However, essentially due to the small de Broglie 
wavelength of the projectile this effect is quite small. These questions will be addressed in Sec. El below. 

2.1 Nonrelativistic (NR) Projectile Velocity 

We discuss briefly the basic ideas of (NR) Coulomb excitation. The condition of no nuclear contact 
is ensured by using bombarding energies below the Coulomb barrier. In the semiclassical theory the 
projectile moves on a classical Rutherford trajectory, see Fig. Ufa), giving rise to a time dependent 
external electromagnetic field. Nonrelativistic semiclassical Coulomb excitation theory is a classical 
textbook example of the application of time-dependent perturbation theory, see, e.g., (201 or EI]- 

The excitation of the target nucleus is due to the static Coulomb interaction, or more generally, (see 
PP, Eq. (II. 1.2)) the mutual electromagnetic interaction between the ions. For NR collisions we can 
neglect retardation, and the electromagnetic interaction can be written as 

W (i, 2) = / / ^ 1 ^/ i(rl)P2(r " 2) ,: /i( ?- /2(r " 2)/c2 , (i) 

J J \T\ — T2\ 

where Pi(fi), ji{rj) are the charge and current densities for the projectile (i = 1) and the target nucleus 
(i = 2), respectively. For non-overlapping charge-densities one can express W(l,2) in terms of the 
electromagnetic multipole moments of the two nuclei. The full expression is given in pp. We give here 



5 



only the most important terms, that is the monopole-monopole and monopole-multipole interactions and 
neglect the less important multipole-multipole interactions. 

We are especially interested in target excitation. In this case we can assume that nucleus 1 is a point 
charge, we have piifi) = Zied(r p (t)) — fi). We can specialize to the coordinate system (r, rp), where 
r is the target coordinate (measured from the center of mass of the target) and r p (t) is the relative 
coordinate between the center of masses of the two nuclei, see Fig. ^a). In the semiclassical theory 
f* p (t) is a classical time-dependent parameter. The charge operator p(r) is given by 

p(f) = eJ2$(r-fl), (2) 

i 

where the r~l denote the proton coordinates. 

Neglecting the multipole-multipole interactions we can write the electromagnetic interaction as 

W(l, 2) = V(l, f p (t)) + V(2, r p (t)) + . (3) 

The last term is the monopole-monopole Coulomb potential between projectile and target which is 
responsible for the Coulomb trajectory of the projectile. The quantity V(2,fp(t)) gives rise to target 
excitations and V(l,fp(t)) to projectile excitations, respectively. 

The monopole-electric-multipole interaction can be obtained from the expansion 

= 471 ^ ^TT^T r v(0^(r p ), (4) 
I' 'p\ Xfi ZA ~r L 'p 

where we have made use of the condition of no nuclear contact so that always r < r p . In this way, the 
variables r and r p are separated. 

The total time- dependent interaction relevant for the target excitations is given by the sum V(2, r p (t)) 
= VE(2,r p ) + VM(2,r p ), where we can write the quantities Ve (electric interaction) and Vm (magnetic 
interaction), see Eqs. (II. 1.12) and (II. 1.13) in [1], as 

V E (2,f p ) = £ i^M(EA,- / i)(-l)V" 1 ^(r P ) (5) 

A>l,/i 

V M (2,f p ) = £ --^{M(MA,- /U )(-l)^r^- 1 LF v (r p ). (6) 

The electromagnetic interaction can be parametrized completely in terms of the electromagnetic mul- 
tipole matrix-elements M(Tr\fi). In the long wave length limit they are defined as 

M(EXfi) = [ d 3 rp(f) r x Y x „(f) (7) 

r x LY^{r) 



M(MXp) = ——- I d 6 rj{r) ■ 



with L = —irx V. Especially we have M(EOO) = Z2e/v47r, where Z 2 is the electric (monopole) charge 
of the target nucleus. 

The condition of no nuclear contact is vital to this approach. In this case we obtain the complete 
separation into the electromagnetic matrixelements of the target. For r > r p the strong interaction 
between the nuclei would come into play and would severely modify the results. This is discussed in 
detail in Sec. El (and also in |2"2"l l2ll]). 
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There are a few dimensionless parameters which characterize the electromagnetic excitation: we 
define the adiabaticity parameter as the ratio between the collision time and the excitation time 

£=— . (9) 

T~exc 

We can estimate the collision time to be t co u = r min /{pfv) and we have r exc = h/AE , where 
AE = hu is the nuclear excitation energy. From this we get 

For NR collisions we have 7 ~ 1, whereas in the relativistic case (see below) the Lorentz parameter 7 
can be much larger than one and the collision time can become very small due to the Lorentz contraction. 
For adiabaticity parameters £ ^> 1 the system can follow the adiabatic ground state and no excitation 
occurs (20] (with probability ~ exp(— £)). This means also that in NR Coulomb excitation one can only 
excite nuclear states for which the long-wavelength limit is valid: Due to the adiabaticity condition 
^ 1 we have ur min <C v. This leads to kR<i <C 1, where k = uj/c and R2 denotes the size of the nucleus 
(see Fig. Ufa)), since r min > R 2 and v < c. On the other hand, for relativistic collisions the collision 
time can be very small and one is able to excite states for which the long wavelength limit is no longer 
valid. 

A second parameter is the strength parameter, which is defined as the strength of the interaction 
potential times its duration (in units of h): 

x = (11) 

Here Vi nt denotes a typical value of the interaction potential. For a multipole interaction of order A, this 
value of the interaction potential is typically e (/|| M(EX) /r m "^; (NR case), the strength parameter 
X is therefore estimated to be 

_ Z ie </||M(£A)||i> 
X ~ five 



rmn 



(The relativistic case will be discussed below: in this case, the collision time has a 7 _1 -dependence; the 
interaction Vi nt is proportional to 7 and % becomes independent of 7.) The strength parameter for the 
monopole-monopole case is the Coulomb parameter, which is given by 

^ 2 e 2 . 

^ = —k • 13 

hv 

In order to obtain the equation describing the time evolution of the system under the influence of 
the electromagnetic field, we expand the wave function of the target nucleus into the complete set of 
eigenstates of the target Hamiltonian $ n with energy hu n . 

(t) = J2 a n(t) exp(-zcj n t)$ n , (14) 

n 

One finds the following system of coupled equations (see, e.g., Eq. (5.29) of [l]) for the excitation 
amplitudes a n (t) 

iha n = ^2 (n\V(t)\m) exp[i(u n - u m )t] a m (t), (15) 

m 

where the time- dependent electromagnetic interaction is denoted by V(t). A formal solution can be 
given using the time ordering operator T as (see Eq. (II. 3. 13) of 0) 

a n (oo) = (n\ T exp f -- J dtV{t)\ |0> , (16) 



7 



where the interaction V is given by 



V(t) = exp(iH t/h)V(t)exp(-iH t/h) 



(17) 



and H is the Hamiltonian of the nuclear system. 

We mention two important limits: if the electromagnetic interaction V(t) is weak enough (strength 
parameter \ <C 1), it can be treated in first order (one photon exchange) and one obtains an expression 
which factorizes into the orbital integrals and the electromagnetic matrixelements. One has pQ 



The orbital integrals S n \^ can be calculated from the kinematics of the process [1] (We have given here 
only the electric part of the excitation, it = E, the magnetic part 7r = M can be found in pQ). Thus 
the experimentally determined excitation probabilities (or cross-sections) are directly proportional to 
the reduced transition probabilities B(n\) (see Eq. (4), p. 93 of [T|). 

Another important limit is the sudden approximation, where the collision time t co u = r min /^v is 
much smaller than the nuclear excitation time r exc = l/oo, i.e., the adiabaticity parameter £ is much 
smaller than one. But the strength parameter \ can have arbitrary values. In this case, we have V ~ V 
and we can neglect the time ordering in Eq. and the electromagnetic interaction can be summed 
to all orders. For large values of £ and x one may use an adiabatic representation of the Schrodinger 
equation, see [T] . It is also often useful to solve Eq. (fT3j) numerically in these cases. (The case of the 
excitation of a harmonic oscillator can be treated exactly analytically, see below.) In contrast to the 
NR case considered in this subsection the velocity for intermediate energy and relativistic collisions v 
is close to c, and the strength parameter x, see Eq. (JT2|) . cannot be too large in these cases and higher 
order effects tend to be small, see also Sec. 12.21 This is in contrast to Coulomb excitation below the 
barrier, where «<Cc and x can be much larger than one. An example is the excitation of a rotational 
band in deformed nuclei pQ, more details will be given in Sec. El 

Although we have shown above that we can treat the electromagnetic interaction in the limit of no 
retardation for NR collisions, one may also treat the electromagnetic interaction in its full relativistic 
form, taking retardation effects into account. This is done, e.g., in [2j using the Coulomb gauge. In 
addition to the contribution from the instantaneous Coulomb interaction (longitudinal photons) one 
gets the one photon exchange contribution due to the H int = j ■ A interaction (transverse photons). 
The multipole expansion of the vector potential A is performed, and it is found 2j that the instan- 
taneous Coulomb interaction is canceled by a certain term (see Eq. f!139|) in the appendix) and only 
the contribution from transverse photon remains. Again there is a factorization into an orbital integral 
and an electromagnetic matrixelement; this time it is the full electromagnetic multipole matrix element 
(see Eqs. (jl41|) and (jl42|0 in the appendix), where the wave number is equal to k = ou/c, i.e., at the 
photon point (please note that the long wavelength approximation is not made in this case). One may 
say that (NR) Coulomb excitation is due to longitudinal (Coulomb) photons, yet it is true that the nu- 
clear structure information is entirely contained in electromagnetic matrixelements at the photon point 
k = uj/c. (For more details see the appendix.) 

There is an important difference of Coulomb excitation to electron inelastic scattering. In inelastic 
electron scattering the photon is virtual (spacelike, i.e., \k\ > ui), and one can study electromagnetic 
matrixelements as a function of momentum transfer. In this way, e.g., the spatial distribution of the 
nuclear charge can be investigated. Also monopole transitions can occur. 

The matrixelements occurring in Coulomb excitation determine the interaction of real (transverse, 
q 2 = 0) photons with the nuclear system. The equivalence of the current and charge matrix-elements 
(see Eqs. (J7J) and ()141J0 in the long wavelength limit is sometimes called "Siegert's theorem", see e.g. 
p. 89ff of [21]. From the above it is clear that the condition of no nuclear contact is the more basic 
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condition than the long wave length limit. In the following discussion of relativistic electromagnetic 
excitation, the long wavelength limit does not generally apply. Yet, due to the condition of no nuclear 
contact, the information about nuclear structure is again entirely contained in the electromagnetic 
matrixelements at the photon point. 

Before going to intermediate energy and relativistic energy electromagnetic excitation let us mention 
the straight line limit of the NR Coulomb excitation theory: From the standard formalism of Coulomb 
excitation using Rutherford trajectories one can also obtain this limit by going to small scattering angles 
9. There remains a characteristic difference between the straight line theory and the straight line limit 
of the Rutherford trajectory: a factor exp(— 7r£/2), see Eq. (II.E.79) of 2J. This factor can be traced 
back to the influence of Coulomb repulsion; the impact parameter b should be replaced by an effective 
impact parameter. This is explained in detail in [2|and |24j . We come back to this in the next section, 
where the relativistic case is discussed. 



2.2 Equivalent Photon Method, Theory of Relativistic Coulomb Excitation 
1. Basic Idea 

The equivalent photon method was introduced by Fermi in a classic paper in 1924 [2HH2H]> see also the 
translation by Gallinaro and White ^ZT\. He considered only the NR case, the relativistic generalisation 
was subsequently given by Weizsacker and Williams [2111211 • This method is lucidly described in the 
textbook of Jackson [Hill, so let us only summarize very briefly the idea: 

We consider the classical straight line motion of a charged particle with charge Ze and impact 
parameter b (see Fig. ^b)). This motion gives rise to a time-dependent electromagnetic field at a given 
point (the target nucleus). The Fourier-transform of the electric field is "compared to the corresponding 
electric field of a suitable spectrum of light" j2Hl- In this way, Fermi introduced the idea of "aquivalente 
Strahlung" (equivalent photons). In this method only the ET-multipole is treated correctly, since the 
spatial variation of the electric field strength over the nuclear volume is neglected. Higher multipoles 
(E2, E3, Ml, M2, ...) show a different behaviour (see, e.g., [7j and below). Higher multipole 
contributions are given correctly only in the limit 7 — > 00 in this approach. 

In the equivalent photon approximation the cross section for an electromagnetic process is written 

as 

r duj 

o — — n{uj)a 7 (u). (19) 
J 00 

where ct 7 (uj) denotes the appropriate cross section for the photo-induced process and n(u) is the 
equivalent photon number. The equivalent photon number in the case of a point particle is given in 
terms of the modified Bessel functions K n as: 



n(u) 



2Z 2 a c 2 

7T V 2 



(20) 



where £ = u}b min /jv. The minimum impact parameter is denoted by b min = Ri + R 2 . Eq. (|2U|) is well 
approximated by 

n(uj) = -Z 2 a\n^- (21) 

7T Ub min 

for v ~ c and 1. 

For £ ^> 1 one can use the asymptotic expression for the modified Bessel functions and obtain a 
spectrum which decreases exponentially with £: 

n(w) » ^exp(-2£). (22) 
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Already from this simple approach the major properties can be read off: the equivalent photon spectrum 
is very soft, it shows essentially an 1/uj dependence. There is an adiabatic cutoff at u — jv/b min . There 
is a (c/f) 2 dependence, which gives a large photon flux for small beam velocities. However, this is 
restricted to small photon energies, u> < jv /b m in- 

2. Semiclassical straight line approximation and exact multipole decomposition 
In the above method equivalent photon numbers of higher (E2, E3, . . . , Ml, M2,. . . ) electromagnetic 
multipoles are in general not treated properly. A beautiful analytic expression which treats all multipoles 
exactly was given in [21], see also |7j. We discuss it briefly now. 

In the relativistic case, retardation can no longer be neglected. One can use, e.g., the Lienard- 
Wiechert potential (see Eq. ([143)) ) in order to describe the (retarded) electromagnetic interaction |24| . 
corresponding to the Lorentz gauge. In the multipole expansion the condition r < r p is used and again 
(see appendix for more details) one can separate the first order amplitude into an orbital part and a 
structure part, which consists of the electromagnetic matrix elements at the photon point (no long wave 
length limit is made). We give only the final formula here: 



a fi = E(- 1 ) m v / 2TTT (~) l G vlm (~) K m (IfMf \ M(ttI - m |J,M,) , (23) 

where the functions G n i m can be expressed in terms of the associated Legendre polynomials, see 

This formula is the basis for the analysis of electromagnetic excitation. The dependence on the 
electromagnetic multipolarity nl, the impact parameter b and the Lorentz factor 7 is clearly exposed. 
We note that the above equation is a first order result, i.e. the strength parameter \ should be much 
smaller than one. In the straight line relativistic case we can write \ as (see Eq. (1.6) of |24J) 

Zl e<f\M(EX)\i> 

x hvb x ' K ' 

We may define an intermediate energy region where jv/c — 7/? = a/7 2 — 1 =~ 1 _ 2. In this region there 
are relativistic effects, see |^TJ , but also effects from the Coulomb repulsion. This region is of special 
interest to us: at GSI we have 7/3 ~ 1-2, at MSU, RIKEN, and GANIL there is £7 < 0.5. 

While in the nonrelativistic theory the Coulomb repulsion of the ions is included in the Rutherford 
trajectory, the problem for the relativistic case (where we have assumed a straight line trajectory) is 
more delicate. A very interesting approach to treat the classical trajectory in a relativistic Coulomb 
problem can be found in [31] [32| . 

It was shown in [21] that a simple redefinition of the impact parameter accounts well for Coulomb 
repulsion: 

b' = b + a -. (25) 

7 

An alternative procedure is suggested in [331 [31] where the following substitution is made: 



b' = a + y/a 2 + b 2 . (26) 

This amounts to replacing the asymptotic impact parameter b by the distance of closest approach 
in small angle scattering. The quantity a® is given in the NR, as well as, in the relativistic case 
by do = hrj/(kv) with hk = jrriaV the momentum of the projectile nucleus. Recently the effects of 
retardation on the electromagnetic excitation in intermediate energy collisions were carefully studied in 



We note some further properties: The high energy 7 — >• 00 limit is investigated in [2l] . The question 
of which value of the magnetic quantum number m gives the largest contribution to a given transition 
depends on the value of k. For small values of k the amplitudes for m = ±1 are largest. The excitation 
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amplitude is proportional to 1/6'. Integrating over b this leads to the logarithmic rise of the cross 
section with 7 only for I = 1, i.e., for dipole transitions (El, Ml). Higher multipoles on the other hand 
tend to a constant in this approximation (i.e. neglecting the m = ±1 contribution). For large impact 
parameters kb > I the largest contribution is due to m = ±1 . This contribution corresponds to the 
Weizsacker- Williams approximation, where the excitation corresponds to the absorption of a photon 
traveling in the beam direction with helicity m = ±1 (see the Eqs. (2.20)-(2.23) in [24J). Now there is 
the logarithmic rise of the cross section for all multipolarities, albeit only for extremely high values of 
7- 

3. Validity of the Semiclassical Method: 

It is well known that for Coulomb parameters 77 ^> 1 the semiclassical approximation can be used. 
This can be seen from the following argument: two different transverse momenta are associated with 
an impact parameter b: The Coulomb repulsion gives a (classical) momentum transfer to the nucleus 
perpendicular to the beam direction. It is given by qr,coui = 2i]h/b. Due to quantal diffraction we 
have a transverse momentum qr,diff = h/b associated with the impact parameter b. The semiclassical 
approximation is appropriate as long as the Coulomb push qr,coui is much larger than the uncertainties 
associated with the wave nature of the particle, i.e., qr,coui 3> qr,diff, or 77 1. A more refined 
derivation of this condition can be found in [SB] and also in App. B of pQ. Another condition for the 
validity of the semiclassical approximation is AE/E <C 1, i.e., the energy loss due to the excitation is 
small. This is well fulfilled due to the adiabaticity condition Eq. (jTUj) . 

From the criterion for the validity of the semiclassical approximation (i.e. 77 ^> 1) we find: for highly 
charged ions with Z\Z 2 > 137 (like the Pb-Pb system) the semiclassical approximation is valid for all 
velocities v, whereas in general we should satisfy the condition Z^Z 2 > 137v/c or v/c < ZiZ 2 /137; for 
sufficiently low velocities the semiclassical approximation becomes always applicable. 

2.3 Quantum Mechanical Treatment of Projectile Motion: PWBA and Glauber Method 

For cases where the Coulomb parameter 77 is smaller than one the semiclassical approximation is not 
valid. Well-known prominent examples are electron or proton scattering on light target nuclei. In these 
cases one uses fully quantal approaches like PWBA or Glauber theory [3112]. The nuclear interaction 
with the nucleus can often be taken into account using the black disk approximation. Characteristic 
diffraction effects arise, see also [371 [TH] and [TH] . 

For 77 ^> 1 one can obtain the semiclassical limit starting from a quantum-mechanical approach 
[T]. In Glauber theory this can be done as follows, e.g., see [38] (see also Eqs. (73)- (80) of [S3])- The 
excitation amplitude has the form (see Eq. (73) of [33]) 



The eikonal phase for a Coulomb potential can be written as x(b) — 2rjln(kb). The inelastic excitation 
amplitude is denoted by a M (cf. Eq. (JTfij) or (fISj)) where \i denotes the ^-component of the angular 
momentum transfer to the nucleus. The integrand in Eq. (J27)) oscillates rapidly as a function of b. In 
the semiclassical approximation we obtain the main contribution for those values of b where the phase 
is stationary. 

We can approximate the Bessel function for large values of qb as 




(27) 



JM h ) ~ ^/2nqb 6XP \ qb ~ ^ 2 ) + 6XP 

We use the stationary phase (saddle point) approximation 




( 



— iqb + iir 



H +1/2 
2 



) 



(28) 
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Figure 2: One-photon excitation e.g. in the (e, e') reaction. In PWBA there is a definite four-momentum 
transfer g M = fc M — k'^ to the nucleus. In contrast to this there is no definite value of q in the Coulomb 
excitation due to the large number of additionally exchanged elastic photons. 



where a phase exp(— iir/4) needs to be chosen for our case, where $"(&o) = ~Q 2 /(^v) < an d which is 
valid for slowly varying functions G(b). The point of stationary phase is denoted by bo, it is determined 
from the condition 

$'(&o)=0. (30) 
$(&) = -qb + 2r]ln(kb). (31) 



With 

the condition Eq. (f3Tf|) leads to 



bo = ^, (32) 
Q 

i.e., the classical connection between impact parameter and scattering angle, or momentum transfer. 
We obtain 

fUfi) « ^J^to)\ eXp ^ (&o) + m ^ a ^ (33) 



Finally we obtain 



da da < 34) 

where P(b ) = 2/ 1 +1 J2 mi ,m f l a /J 2 ; analogous to the expression in the nonrelativistic case [21 II]- 

The main difference of electromagnetic excitation in hadron-hadron scattering, as opposed to electron- 
(lepton-) hadron scattering, is due to the strong absorption at impact parameters less than the sum of 
the two nuclear radii. Thus, for 77 ^> 1 the semiclassical method is very suitable to take these effects 
into account, because the impact parameter is the relevant variable and not the momentum transfer. 
Neglecting strong absorption effects it can be shown that (first order) SCA and PWBA give identical 
results for total cross sections (where the limit of the momentum transfer q ma x must be extended to 
00; this is in general a good approximation as the contribution to the cross section for large q is rather 
small ). For the nonrelativistic case this is shown, e.g., in [SU], the relativistic generalization is given by 
[7j. Likewise the equivalence of the coupled channel semiclassical theory and Glauber theory (i.e., with 
the electromagnetic interaction treated to all orders) is given in |3U]. It should be noted however that 
differential cross sections are different in the PWBA and SCA approaches, and thus the appropriate 
theory should be used in each individual case. 

In the PWBA there is a definite four-momentum transfer q^ = — k' , that is, a definite energy 
transfer uj = q and a definite three-momentum transfer q = (q~r,qi). To a good approximation qi is 
given by the minimum momentum transfer qi = q m i n — u/v for small angle scattering and small energy 
loss. The invariant mass of the photon is Q 2 = —q 2 = q\ + (uj/('~fv)) 2 . We always have Q 2 > 0: the 
exchanged photon is virtual (spacelike). 
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In contrast to this, the photon, which is exchanged in Coulomb excitation, contains a sum over 
virtual photon momenta; they conspire in such a way that only an electromagnetic matrix-element 
survives, which corresponds to the interaction with a real photon Q 2 = (see appendix for details). 

In principle, heavy ion electromagnetic excitation can also be treated within field theory. This is 
however rather complicated and has to our knowledge never fully been done in practice. Due to the 
validity of the semiclassical approximation it is also not necessary to do so. Let us indicate how one 
would proceed ^TJ|321 : There are infinitely many graphs which have elastic photon exchanges between 
the two ions. For small angle scattering and high energies one can make use of an relativistic eikonal 
approach [33J 133]. In this approach infinitely many Feynman graphs are summed. The result can 
again be expressed in terms of an expression in impact parameter space with a relativistic eikonal 
X- The semiclassical approach can be recovered in exactly the same way as already explained above. 
From the saddle point (stationary phase) approximation one obtains again the relation b = 2r]/q. For 
b < b m i n = R1 + R2 there is nuclear contact and and graphs which involve the strong interaction have to 
be added, see jH]. In the case of nucleus-nucleus collisions the black disk model can be applied. This 
means that the eikonal x is a very large imaginary number for b < b m i n = R\ + i?2- 

3 Higher Order Effects 

The importance of higher order effects is governed by the strength parameter x- This parameter is 
proportional to 1/v, see Eq. (II 2|) and Eq. (|24)) above. For slow collisions the Coulomb field acts for a 
long time and this parameter can be quite large. E.g., rotational bands in deformed nuclei can be very 
well excited in this way and this has been a fruitful field of research PQ. For collisions with v ^ c higher 
order effects tend to be small. The strength parameter x will soon reach its lower limit, which we have 
for v = c in Eq. JEJ) and Eq. (1231) . 

Let us deal in the first subsection with a very interesting higher order electromagnetic excitation 
process at intermediate and relativistic energies: the excitation of the Multiphonon Giant Dipole Res- 
onance, especially the Double Phonon Giant Dipole resonance (DGDR). In the following subsection we 
will discuss a more disturbing effect: corrections of one-photon exchange effects due to multiphoton 
exchange (a similar problem exists for (e, e') scattering named dispersion corrections see, e.g. |21j). 

3.1 Multiphonon Giant Dipole Resonances 

The Giant Dipole Resonance (GDR) is a collective mode par excellence. This can directly be seen from 
the fact that the energy-weighted TRK-sum rule (Eq. ()1U9|) below) is exhausted by this mode. With 
its large B(El)-value this transition is very strongly excited electromagnetically in intermediate energy 
heavy ion collisions. For heavy systems the excitation probability in grazing collisions is of the order of 
30 percent. From this one can conclude directly that higher order effects must play a role. This was the 
situation in the middle of the eighties, and the question was to which states does the GDR couple most 
strongly. There is the Axel-Brink hypothesis that a GDR is built on any nuclear (excited) state. In this 
spirit a multiphonon harmonic oscillator model was adopted in [T21 I45] and higher order electromagnetic 
excitation was calculated within this model. In the harmonic oscillator model the interaction can be 
included to all orders and a Poisson distribution of the multiphonon states is found, as is well known, 
see, e.g., [T]. We give some key steps in the derivation of this result, see also [46]. 

In terms of the corresponding creation and destruction operators at and a the Hamiltonian of a 
harmonic oscillator is given by 

H = hu){(fa+-), (35) 
where u denotes the energy of the oscillator. We have the boson commutation rules [a, at] = 1 and 
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Figure 3: One- and two-photon excitation of the DGDR. The one-photon excitation to the DGDR is 
forbidden (0 + — > + ) or strongly hindered (0 + — > 2 + ). 



a, a 



0, respectively. Only one mode is shown explicitly, in general one has to sum (integrate) 



over all the possible modes. 

We assume that the interaction V is linear in the creation and destruction operators: 

V{t) = f{t){a + at). 



(36) 



In this case we can calculate V (Eq. (|17j)) explicitly using the boson commutation rules given above 
and the expansion 



One finds 



V(t) = V(t) + it[H , V(t)} + ^-[H , [H , V(t)}} + 



V(t) = f(t) (aV^* + ae lu)t 



(37) 



(3* 



Now we can convince ourselves that the commutator of V at different times t and t' is a pure c-number. 
In this case we can disregard the time ordering operator in Eq. (J16J) and obtain an exact analytical 
answer, up to an unimportant overall phase factor (for details see, e.g., [T]). 

This leads to the excitation of a so-called coherent state, see |47| UK] . For the excitation of multi- 
phonon states this is explicitly shown, e.g., in [JH]- One has 



n i-i(«*at+wa)|g 



-iu*r _i, 

-e 2' 



where u is the c-number 



u 



dtf(t) exp(-iujt) 



(39) 



(40) 



In order to show this the operator identity e A+B = e A e B e~^ A ' B ^ was used, which is valid for two 
operators A (= — iu*(v) and B (= —iua) for which the commutator is a c-number. 
The excitation probability Pat(6) of an iV-phonon state is given by 



P N (b) 



P( 1 )(6) 



JV 



exp 

IvT 



-pW(6)" 



(41) 
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Figure 4: (upper panel) Integrated experimental cross sections for the one-phonon GDR and the DGDR 
in 208 Pb obtained from different targets of nuclear charge Z. In case of the DGDR the calculated values 
are multiplied by a factor of 1.33 to fit to the experimental data, (lower panel): ratio of the experimental 
cross sections for the DGDR in 208 Pb to the harmonic approximation. The mean value and its error are 
indicated by solid and dashed lines , respectively. Taken from [HU], Fig. 2 there, where further details 
can be found. 



where P^\b) is the excitation probability of the GDR for an impact parameter b calculated in first 
order perturbation theory. It is given by [7j 



P«(6) 



2a 2 ZfN 2 Z 2 1 



(42) 



where the energy uj of the resonance is given by uj = 80A _1 / 3 MeV and where A^ 2 and A 2 are the neutron 
and mass number of the nucleus, which is excited. 

Certainly the most difficult question at that time was the width of the multi-phonon states (It still 
is today.) 

Subsequently detailed experiments [HUE2] were carried out at GSI and the two-photon excitation of 
the DGDR was clearly observed in a two step excitation process, see Fig. 01 Generally, the experiments 
also confirmed the simple theoretical approach. In Fig. 01 we show the dependence of GDR and DGDR 
excitation on the charge number Z of the target. For the GDR we have a one-photon exchange, i.e., 
there is a Z 2 dependence. For the DGDR we have a two-photon exchange, i.e., a Z 4 dependence, 
essentially in agreement with the experimental result |oT?j . 
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However, interesting differences to the harmonic model showed up, which motivated various groups 
to do more detailed theoretical calculations. Most conspicuously, the width of the DGDR tends to be 
somewhat less than twice the width of the one-phonon GDR. The DGDR is a remarkable example of a 
large amplitude collective motion with damping. We refer in this respect to the reviews jTT] and |53j . 
where also further references can be found. The influence of the damping on the GDR- and DGDR- 
excitation cross sections was investigated in |54j . a more sophisticated approach was adopted in j55j . 
This work also puts the approach of Carlson et al. [SHI EZ| on a firmer theoretical basis. A schematic 
model which treats the coherent (collective) excitation modes and their coupling to more complicated 
states was considered in 

It would also be very interesting to study the DGDR in deformed nuclei like 238 U. What will the 
splitting of the DGDR be due to the nuclear deformation? What is the influence of the nuclear shape 
on the minimum impact parameters? Clearly, a closer look at these problems would be very interesting 
in the future. 

Due to the very collective nature of the GDR, the electromagnetic couplings in this case are very 
strong. As we have seen in this section, interesting effects follow from this. Even at the relativistic 
heavy ion colliders RHIC and LHC the GDR, and possibly also DGDR, play an interesting role, which 
is briefly discussed in Sec. 16.21 

Finally we mention that multiphonon states were also reviewed in jSHl, where emphasis is also put on 
the nuclear excitation mechanism (which becomes important for lower beam energies) and quadrupole 
phonons. Multiphonon states were first investigated in pion-charge exchange reactions [60j . 



3.2 Corrections to One-Photon Exchange 

Now we turn to cases where the higher order effects are less pronounced. But still they are there and 
they are essentially a correction to the most interesting dominant first order effect. There are various 
methods to take higher order electromagnetic effects into account theoretically. If the coupling is strong, 
a coupled channels approach is appropriate (see Eq. (|15Jl ). sometimes higher (especially second) order 
perturbation theory is sufficient. There is a sum over all intermediate states n, which are (considered to 
be) important pQ. Another approach is to integrate the time-dependent Schrodinger equation directly 
for a given model Hamiltonian [SH E21 EH1 IMj • In this way all orders are taken into account. 

If the collision is sudden (£ <C 1), one can neglect the time ordering T in the usual perturbation 
approach, see Eq. (|16j) (sudden approximation) . The interaction is again summed to infinite order. In 
order to obtain the excitation cross section one has to calculate the matrix element of this operator 
between the initial and final state. Intermediate states n do not appear explicitly. 

In order to see the typical effects and how they depend on the relevant parameters we now study 
higher order effects in the simple model of ^3]. Let us briefly describe this model. It will be discussed 
also in the following section, where this model is discussed in the so-called post-form DWBA (distorted 
wave Born approximation). (The method used here is more closely related to the prior-form DWBA.) 

We assume a target nucleus with charge Z (to avoid unnecessary complications we assume an infinite 
mass for a pedagogical discussion, this is however not essential). We want to describe the breakup of 
a bound state a = (c + n) with binding energy Eq in the Coulomb field of this target. There is 
the Coulomb interaction V c = Z c Ze 2 /r c between the target and the core c. We assume a zero range 
interaction between the core and the neutron. This can be viewed as a deep square well potential with 
depth V extending up to a radius a with a 2 V held constant when a — > 0. In this potential an s-wave 
bound state can exist. 

In order to calculate the breakup process Z + (c + n) — > Z + c + n we assume that the (c + n)- 
system moves on a straight line trajectory with velocity v and impact parameter b in a semiclassical 
model. We restrict ourselves to electric dipole transitions. The El and E2 effective charges are given 



by Zty f = Z c I = 1, 2, • • Due to the smallness of Z£> f it is a good approximation to neglect 



'eff 



(2) 



eff 
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the E2 transitions for the neutron-core case (This is in contrast to the proton-core case to be discussed 
below). Analytical results were obtained for I s * and 2 nd order electromagnetic excitation for small 
values of the adiabaticity parameter £. We are especially interested in collisions with small impact 
parameters, where higher order effects tend to be larger than for the very distant ones. In this case the 
adiabaticity parameter £ is small. For £ = (sudden approximation) we have a closed form solution. 
(This is a special example of the general formula, see Eq. (fTo^) but without the time-ordering operator). 
In Eq. (37) of pH] the angle integrated breakup probability is given. We expand this expression in a 
scaling parameter 

2ZZ c e 2 m n m n ri . . 

nv(m n + m c )bK [m n + m c )bK 

where the parameter k is related to the binding energy Eq by 



and the reduced mass /i is given by 



In leading order (LO) we obtain 



E ° = w (44) 

A* = (45) 
m n + m c 

This parameter is directly related to the strength parameter x- We define another scaling param- 
eter 

x = ^, (46) 

K 

where the wave number q re i is related to the energy E re i of the continuum final state by 

E Tcl = ^fsL. (47) 

= igy ^ (48) 
(ig re i 37tk (1 + x 2 ) 4 

The next to leading order (NLO) expression is proportional to y A and contains a piece from the 2 nd 
order El amplitude and a piece from the interference of 1 st and 3 rd order. We find [65J 

dP NLO _ _16_ 4 x 2 (5 - 55x 2 + 28x 4 ) 

dqrel ~ 37TK y 15(1 + X 2 ) 6 ' ^ ^ 

These are very simple and transparent results for the Coulomb dissociation of neutron-halo nuclei. 
There are two scaling parameters x and y. The parameter x controls the shape of the relative energy 
distribution, whereas the parameter y plays the role of the strength parameter. We have given these 
formulae explicitly, since they show most of the important features in a simple analytical formula: as is 
to be expected, y is proportional to the target charge Z and the effective El-charge Z c m n / (m n + m c ). It 
is proportional to 1/v, again this is the special case of a general result. We also have the 1/6 dependence 
characteristic of El transitions. 

The integration over x and the impact parameter b can also be performed analytically in good 
approximation, for details see [65] . We can insert the corresponding values for the Coulomb dissociation 
experiments on n Be and 19 C [HZ1 EHJ in the present formulae. We find that the ratio of the NLO 
contribution to the LO contribution in the case of Coulomb dissociation on 19 C jBE| is given by —2%. 
This is to be compared to the results of jHHI, where a value of about —35% was found (A further 
discussion of this is found in [65J). A version of the theoretical model used in [HE], see a ls° is 
discussed below in Sec. |U 
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Figure 5: Differential cross sections integrated over scattering angle from 0° to 3° for the Coulomb 
dissociation of 67 AMeV 19 C scattered on 208 Pb as a function of the relative energy. The solid line 
shows the results of the LO expression (Eq. (|48p ) . the dashed one includes also the NLO contribution 
(Eq. (|49))). The dotted line shows the results of a semiclassical calculation for El only in first order, 
the dashed-dotted line (which almost coincides with the (dashed) LO+ NLO results) for the numerical 
solution of the Schodinger equation with El. Both I s * order calculations (solid and dotted line) are 
also in close agreement with each other. The points are the experimental results from [HHl- Reproduced 
from Fig. 3 of [S3]. Copyright (2001) by the American Physical Society. 



In Fig. we show a comparison of theoretical and experimental differential cross sections for Coulomb 
dissociation of 19 C , see Fig. 3 of jHS]- In this figure, LO and NLO results (Eqs. (]4"%j) and (JUIJ)) are also 
compared to full dynamical calculations (solution of the time-dependent Schrodinger equation). 

Let us now deal with the breakup of a proton-core system as another typical case. Examples 
are 8 B = 7 Be + p or 17 F = 16 + p. In this case, the .El-effective charge is small and the in- 
effective charge is large, as compared to the neutron-core case. The effective charges are given by 
Zfjj = Z c { y ^ rr ^ v m ) l + Z p ( m rn ^ m ) l (I = 1,2,...) where m p and m c are the masses of the proton and 
the core, Z c is the charge number of the core and Z p = 1. Also in view of the importance of the 
8 B Coulomb dissociation for the solar neutrino problem (more see below in Sec. |7J this case has been 
studied extensively with all the methods available, which we now briefly survey. In most models, some 
nuclear structure quantities enter and it is a pity that somewhat different assumptions (for example 
about El or E2 matrix elements) are used in the different calculations by the various groups. This 
makes the comparison of the different reaction calculations somewhat more difficult. How does one 
separate the effects from the nuclear structure input and the effects of the different approximations in 
the reaction theory? Yet, quite a conclusive picture has emerged. Using perturbation theory, higher 
order effects in the Coulomb dissociation of 8 B were studied in |14[ I15j. This approach works best for 
the high beam energies. Results are given in Fig. ESI see Sec. UJ below. 

The CDCC (Coupled Discretized Continuum Channels) method was also applied, for recent refer- 
ences see |70| ITT] . The effects of E2 and higher order excitations are included. Since the relative motion 
of the projectile and target, with its huge number of partial waves is treated in a quantum mechanical 
way, this kind of method becomes numerically more and more involved for high energies. In it was 
seen that (not unexpectedly, since the corresponding 77-value is much larger than 1) the semiclassical 
approximation compares well to corresponding fully quantal calculations at sub-barrier energies. This 
is a great simplification in the numerical evaluation. At higher energies, where also nuclear (grazing) 
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Figure 6: Diffraction dissociation probabilities of the F ground state at 10 and 40 AMeV on a 208 Pb 
target. The points are the results of the full dynamic calculation. The long-dashed and steep solid line 
is the result of an eikonal calculation with and without core-absorption effects. The flatter solid line is 
the result of a first order calculation including El and E2 contribution. The short dashed line is the 
result of the dynamic calculation including only the Coulomb interaction. Taken from [22], Fig. 5 there. 



collisions can become important, a Glauber approach is very useful |3I]. 

Instead of expanding the wave function into the nuclear basis states one may also study the time- 
development of the nuclear wave function directly. With modern powerful computational methods it is 
now possible to solve the time-dependent Schrodinger equation very efficiently [73J 1221 E2 EH E3 • It 
may be considered to be a drawback of these approaches that they rely on choosing a specific (simple) 
nuclear Hamiltonian, like a particle-core model with a phenomenological interaction between them. 
However, it will cover the essential aspects of the problem in most cases. In practice such approaches 
turn out to be very useful and also give us some general insight into, e.g., the importance of higher 
order effects. 

Esbensen and Bertsch j22] investigated the two-body (p+ 16 O) breakup of 17 F on Ni and Pb targets 
at 10-40 AMeV. Higher order, as well as, nuclear effects are important and we report on some of their 
results now. A "dynamic polarization (Barkas) effect" is found. It is essentially due to the interference 
of a (first order) El amplitude and a (second order) E1-E2 amplitude. The effect is proportional to 
Z 3 . In Fig. IHlwe show the breakup probabilities of 17 F on a Pb target at 10 and 40 AMeV as a function 
of the impact parameter. 

In P one may find various approaches to study deviations from the sudden approach, see, e.g., the 
discussion in Sec. II. 3. It would be of interest to further explore this in the future. 

In summary we may say that higher order effects are well under control. After all we have a fair 
knowledge about nuclear structure and a very good knowledge of how to treat the electromagnetic 
interaction. Due to present days' computing power many quite involved calculations have become 
possible and they are essential for the analysis of the experimental data. Simple models and analytical 
solutions are also useful: they show the dependence on the various parameters. 
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4 Analytically Solvable Model for the Coulomb Breakup of 
Neutron Halo Nuclei 



Breakup processes in nucleus-nucleus collisions are complicated, in whatever way they are treated. They 
constitute at least a three-body problem, which is further complicated due to the long range Coulomb 
force. Exact treatments (like the Faddeev-approach) are therefore prohibitively cumbersome. On the 
other hand, many approximate schemes have been developed in the field of direct nuclear reactions, and 
these approaches have been used with considerable success [THj. In this context we wish to investigate 
a realistic model for the Coulomb breakup of a neutron halo nucleus. This model has already been 
mentioned in the previous section. With the operation of exotic beam facilities all over the world, these 
reactions (previously restricted essentially to deuteron induced reactions) have come into focus again. 
The Coulomb breakup of these nuclei is of interest also for nuclear astrophysics, since the breakup cross 
section can be related to the photo-dissociation cross section and to radiative capture reactions relevant 
for nuclear astrophysics jTH], see Sec. below. 

In this section we study the Coulomb breakup of a neutron-halo nucleus with a zero-range neutron- 
core interaction. A full quantum mechanical treatment can be done with two different treatments of 
the final state. In the so-called prior-form the Hamiltonian H is split in the initial and final channel 
into 

H = H p / ior + Vf lor , (50) 

where 

Hf or = T + V a + V on (51) 

and 

yprior = y c + y n _ (52) 

where Vi denotes the interaction of % = c, n with the target nucleus and V cn is the relative core-neutron 
interaction. The interaction V a is an optical model interaction of the a = (c + n) bound state with 
the target nucleus. The final state is a product of the wave function of the center of mass (with the 
coordinate vector R) and the relative wave function of the unbound (c + n) system [ZE] • 
In the so-called post-form the Hamiltonian H is split in the final channel into 

H = Hf st + Vf os \ (53) 

where 

H f st = T + V n + V c (54) 

and 

Vf ost = V nc . (55) 

The interactions V n and V c are treated to all orders in this way, whereas the relative interaction between 
c and n is only treated in first order. In the previous section we have studied the prior-form DWBA 
(in its semiclassical limit, i.e. for large values of the Coulomb parameter 77). Accordingly we have used 
a decomposition into a cm. motion of the projectile system and its relative wave function in the initial 
and final channel. 

An important benefit of the present model is that it can be solved analytically in the CWBA 
(neglecting the nuclear interaction of the projectile with the target we can call the DWBA the Coulomb 
Wave Born Approximation CWBA), see Sec. 14.21 Analytic expressions also exist for the PWBA limits 
of both prior- and post-form, which we can compare with each other. Thus this model constitutes 
an ideal theoretical laboratory to investigate the physics of breakup reactions. We can study certain 
limiting cases and its relation to the semiclassical approximation, which is mainly used in the analysis 
of experiments. Especially the effect of postacceleration (to be explained in more detail below) can be 
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studied in a unique way in this approach. We will first study the plane wave limits of both the prior- 
and post-form and compare them with the semiclassical limit. The question of higher-order effects will 
then be discussed in the full CWBA approach. 



4.1 Theoretical Model, Plane Wave Approximations and Qualitative Discussion 

We consider the breakup of a particle a = (c + n) (deuteron, neutron-halo nucleus) consisting of a 
loosely bound neutral particle n and the core c (with charge Z a = Z c ) in the Coulomb field of a target 
nucleus with charge Z 

a + Z -> c + n + Z. (56) 

To simplify the kinematical relations we assume in this section that the target is infinitely heavy (this 
is not necessary, see below). We assume that the a = (c + n) system is bound by a zero range force, 
see Sec. EJ The potential V cn is adjusted to give one s-wave bound state with a binding energy E . 
Neglecting the nuclear interaction of c and n with the target {pure Coulomb case) the Hamiltonian of 
the system is given by 

ZZ e 2 

H = T n + T c + — + V m (r). (57) 

' c 

The bound-state wave function of the system is given by 

7Texp(- K r) 



where the quantity k is related to the binding energy E , see Eq. (j44j) above. 

The incoming particle a has a momentum q a and the momenta of the outgoing particles c and n are 
denoted by q c and q n , respectively; see also Fig. [7| below. It is useful to introduce also the relative and 
cm. momenta q re i and q cm which are related to q c and q n in the following way: 

Qrel = ; (59) 

m n + m c 

q C m = q c + q n - (60) 
These momenta are further constrained by energy conservation 

?a K _ qc _|_ Qn _ Qcm |_ ffrel (61) 



2m a 2 ix 2m c 2m n 2(m c + m n ) 2ii 

From the experimental point of view it is interesting to mention at this point a magnifying glass 
effect for a given scattering angles of the particle c and n there is a minimum relative energy E™™. 
The relative energy E re i = q^ el /(2Li) is a small quantity, close to E™™. It can be determined accurately 
from the measurement of the momenta q c and q n (which are both large), see, e.g., JE] and Fig. 4 of 

It is interesting to write down the amplitude for the breakup reaction in Born (plane wave) approx- 
imation. This approximation is valid for small values of the Coulomb parameter rj a . This parameter 
characterizes the strength of the Coulomb interaction. In the applications considered here r\ a is usually 
much larger than one, but still, this example will expose many characteristic features of the breakup 
process. In the prior-form we split the Hamiltonian, see above, into H = H Q + V c with H = T + V cn 
(i.e. we we set V a = ). Accordingly we have in the initial state a plane wave for the cm. motion and 
a bound state wave function; in the final state we have again a plane wave of the cm. motion and the 
relative continuum wave function, which in the zero-range approximation is given by 

\q re i) = exp(iq rel ■ r) - - — — — , (62) 

(k + iq re i)r 
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i.e., only the s-wave is modified by the short range potential V cn . Using the Bethe integral 



1 1 1 r ex P 



R-^r 



iq- (R— —r 



d A q 1 \ ma \ (63) 



2vr 2 J * q 2 



we can separate the matrix element into the relative and cm. coordinates. The T-matrix in the Born 
approximation (PWBA) is found to be 

77 p 2 

T Born = A7r £^_ (A } (64) 

q z 

where Ap is related to the momentum transfer (or Coulomb push) to the target nucleus 

Q = Qa~ Qcm = (fa- Qe ~ Qn (65) 

by Ap = (m n /m a )q. This Coulomb push has a perpendicular component q± and a component in the 
beam direction qn. For high energies we have qn = uo/v (corresponding to the minimum momentum 
transfer). The amplitude a/j can be calculated analytically, see, e.g., Eqs. (33), (34) of [13 . It is found 
to be 

a f i = V8iTK,(a F T + as)- (66) 
The quantity oft is essentially the Fourier transform of the Yukawa wave function, given by 

1 

oft = (67) 

(q rel - Ap) 2 + k 2 

and as takes the s-wave scattering part of the continuum wave into account: 

i(K + iq re i) k + i(q re i + \ Ap\ 

as = — m — t~i — ™ ?tt In ; h — rr- 68 

2\A P \(k 2 + q 2 rel ) K + t(q rel -\Ap\) 1 ' 

We are interested in forward angle scattering at high energies. The value of qn is usually small and 
we have q ~ q a 6 where 9(<g. 1) is the scattering angle. We can introduce the modulus of the elastic 
scattering Coulomb amplitude 

/- = ^ » % (69) 
q 2 q a v 2 

We find T = (27r/m a )/ cou ^ajj. For very small values of q, or equivalently Ap, the terms a F r and a s 
nearly cancel each other (the bound and the continuum wave functions are orthogonal) and the dipole 
term will be dominant. We find 

a fi (q reh Ap) w V&TK 2 f el A ^ (70) 

We can compare this expression with two others: the semiclassical straight line expression jT^] in Sec. H3 
above and the Born limit of the post form CWBA. 

First, we recall the results of [13] and Sec. El above. In this approach the semiclassical straight 
line approximation is used. It is valid for high beam energies and for rj a ^> 1. In the semiclassical 
approach (corresponding to the prior form decomposition of H), one now calculates (impact parameter 
dependent) breakup amplitudes a/i(6); the breakup amplitude is given by 

f breakup fcouiafi (71) 
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and the cross section is 

dcr/dtt = da ruth /dQ \a fi (b)\ 2 , (72) 

where 

\fcoui\ ■ (73) 



dQ 

The impact parameter b is related to the Coulomb push by the semiclassical relation b = 2r] a /q. We see 
that the formula for the Born approximation is the same as the one derived for the semiclassical sudden 
limit. Please keep in mind however that the ranges of validity of the two approaches, do not overlap: 
for the Born approximation we have rj a <C 1 while the semiclassical approximation requires i] a ^> 1. In 
the sudden limit we have to = and there is only a transverse momentum transfer q±. The breakup 
amplitude in the sudden limit is given by, see Eq. (33) of [T3] 

a fi = (q\exp(iKp- r) |0) . (74) 

It should be kept in mind that in ^3] only the electric dipole component is considered. The above 
equation leads to 

a fi = i&p ■ (q\ f |0) (75) 

in lowest order of Ap or q. This expression is the same as Eq. (ffH|) . In this case there is no interaction 
between the c and n in the final state p— wave. In higher orders of q, however, this interaction will show 
up. 

Now we compare these prior form expressions to the post-form DWBA scheme. In the prior form we 
have a c.m.-motion (either classical or quanta!) of the center of mass (with coordinate vector R) of the 
(c + n)-system (either in the ground- or excited state) in the Coulomb field ZZ c e 2 / R. If this potential is 
neglected because of the high projectile energies, one has a plane wave or a straight line trajectory. The 
relative wave function of the (c — n)-system is governed by the internal Hamiltonian H cn = T cn + V cn . 

In the post-form the Hamiltonian in the final state is split as in Eqs. (|53|) - (|53j) . I.e., we take the 
Coulomb interaction between Z and c fully into account in the final state (i.e., postacceleration effects 
are fully considered) and treat V cn to lowest order only. In our model there is no resonance structure 
in the c + n continuum. We can expect that V cn is not so important. This is clearly a good assumption 
for the deuteron and will also hold for other neutron halo systems. Similarly a cm. Coulomb wave 
function is used for the initial state. In the final state the broken up system is treated in a different 
way, more appropriate if there are strong postacceleration effects. This is, e.g., found to be the case for 
low energy deuteron breakup, see [77J EE]- For this see the next subsection. 

We first give explicitly the (Born) plane wave limit |7H] of the post-form DWBA (or CWBA). The 
Born approximation in this case is obtained by substituting for the full Coulomb wave functions in 
Eq. ()85|) below their I s * order expansions in V c . We find (see Eq. (3) of [7JJ]) 



T = D f cou i4ir 



1 m r 
+ 



Ql - (<fn + <?c) 2 m a {ql - (q a - q n ) 2 ) 
where the zero-range constant is given by 



(76) 



h 2 



D = ^V8^. (77) 

As above, f cou i is the modulus of the Coulomb scattering amplitude Eq. (jHHj) . 

This expression shows a similarity to the Bethe-Heitler formula for bremsstrahlung Z + — ► 
Z + e~ + 7. The Bethe-Heitler formula has two terms which correspond to a Coulomb interaction of the 
electron and the target followed by the photon emission and another one, where the photon is emitted 
first and then the electron scatters from the nucleus. Now we have a Coulomb scattering of the incoming 
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Figure 7: The two bremsstrahlung type of graphs, which describe the Coulomb breakup in the post- 
form Born approximation. Three-momentum conservation at each vertex determines the intermediate 
momenta q^. 

particle followed by breakup a — (c + n) — > c + n and another term, where the projectile a breaks up 

into c + n, and subsequently, c is scattered on the target Z. 

The two terms in the parenthesis correspond to the two graphs shown in Fig. [7| 

It is useful to rewrite the denominators of Eq. ()76|) in terms of the relative and cm. variables q re i 

and Ap (Eqs. JBDJ) and flUJ)) as 

^-d, = ^M^ + * 2 ) (78) 

m n m c 

and 

ql - & - q n ? = -I^±^V + (Qrei - Ap) 2 )), (79) 
where we have used the energy conservation Eq. (J61)) to obtain this result. For the T-matrix we obtain 

T = y/^M^L( 1 i ). (80) 

m a "qli + K 2 (q rel -Ap) 2 + K^ 

For small values of Ap, or equivalently q, the two terms almost cancel. 

We can expand this equation in the small parameter Ap. In lowest order we obtain 

2tt 2Ap ■ q rei 



T = V&^f mll - WT ^. (81) 

We can compare this to the limit of small Ap in the prior form PWBA, Eqs and f)70j) . We see that 
these limits agree with each other. 

For larger values of q characteristic differences will appear: in Eq. the n-p final state interaction 
in the / = channel is treated properly. This is not the case for Eq. ()76|) or (jHUj) . In this equation, on 
the other hand the Coulomb force on the core c in the final state is treated properly. 

In the Born approximation it is only the momentum transfer q '= q a — q c — qn and not the momentum 
q a itself which enters into the matrix element. In the semiclassical approximation, on the other hand, the 
quantity v a = q a /m a appears explicitly. This quantity enters in the adiabaticity parameter £, a concept 
not present in the plane wave approach. The (nonrelativistic) adiabaticity parameter (see Eq. (fTU|) ) is 
given by 

2rj a u/v a 2-q a q\\ 

£ = Lob/v a = = ^ (82) 

Q q± 

We have q\\ = u/v a and, for not too large impact parameters, q±_ ~ q. As we know, for large values of 
the adiabaticity parameter, i.e. when 2rj a q\\ 3> q±, the breakup probability is suppressed exponentially 
with £. In the next subsection we will also encounter a simple limit of the full CWBA expression, very 
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much related to the plane wave (Born) treatment of the post form CWBA given above. It is valid for 
small values of £ and arbitrary values of rj a . We will also treat the case of large values of r\ a and £ there. 

It is an interesting and difficult task to specify more clearly the ranges of validity of the post and 
prior formulations. In the simplest high energy limit we have just found that there are regions where 
both approaches agree with each other. For low energy deuteron breakup the post-form CWBA is very 
reasonable, whereas the prior form is very poor jHOj- What about high energies? In second order of Ap 
differences show up. Which of the two approximations is reasonable, and in which region of parameter 
space? The relevant parameters are x, y, and £, what is the range of validity of x- and ^/-scaling? In 
this context we mention a recent paper [HI], where final state three-body Coulomb effects are discussed 
for the 208 Pb( 8 B, 7 Be + p) 208 Pb Coulomb breakup reaction. 



4.2 CWBA 

In the post-form CWBA the Hamiltonian is split in two ways corresponding to the initial and final 
states. We have H = Hi + Vi where 

ZZ e 2 

Hi = T + V nc (r) + V a (R) = T R + T r + V nc (r) + (83) 

i.e., the initial state separates in the r and R coordinates. In the final state we have H = Hf + V) with 

ZZ e 2 

H f = T + V c (r c ) = T rc + T Rn _ {Ac) + — c —. (84) 

' c 

The T-matrix for the reaction Eq. 1)56)1 can be written as [H2j 

T = (xi ] ^ n \Vnc |x£Vo) = D J d'Rx^^e-^x^iR), (85) 

with the zero range constant D as given in Eq. (|77)l above. 

In order to compare with experimental data, one has to take the finite target mass wia into account 
in a standard way. For this one replaces the momenta of the previous section (where the target mass 
was assumed to be infinite) by the ones with the correct kinematics and the masses m a , m c , and 
m n would have to be replaced by the reduced masses (m a m a) / f {m a + ttia), {jn c rnA){yriA + rn c ) and 
{m n {mA + m c ))(m n + m c + m^), respectively. Note that if we would have included also an interaction 
V n such a separation would not have been possible. 

The initial state is given by the incoming Coulomb wave function with momentum q a and the 
halo wave function O . The final state is given by the independent motion of the core described by the 
outgoing Coulomb wave function xk^ m the Coulomb field of the target nucleus Z with asymptotic 
momentum q c and the free neutron with momentum q n , described by a plane wave. 

In these wave functions the Coulomb interaction is taken into account correctly to all orders. The 
present post-form description, Eqs. ()85)) . includes therefore the effects of postacceleration. Postacceler- 
ation can be viewed as a higher order electromagnetic effect, see Sec. EH and refers to the fact that (at 
low beam energies) the core c has a larger final state energy than what one would get from sharing the 
kinetic energy among the fragments according to their mass ratio. Postacceleration arises in a purely 
classical picture of the breakup process. This is nicely discussed in |53] (We show their Fig. 5 here as 
our Fig. IHJ). The nucleus a = (c + n) moves up the Coulomb potential, loosing the appropriate amount 
of kinetic energy. At the breakup point (marked as "breakup occurs here", see Fig|HJ), this kinetic energy 
(minus the binding energy) is supposed to be shared among the fragments according to their mass ratio 
(assuming that the velocities of c and n are equal). Running down the Coulomb barrier, the charged 
particle c alone (and not the neutron) gains back the Coulomb energy, resulting in its postacceleration. 
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a=(c+n) 
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c 



Figure 8: A schematic classical view of Coulomb-breakup, as adapted from jH3]. The distance from the 
target nucleus to the breakup point is denoted by r. 



Figure 9: Comparison of calculations and measurement for the deuteron breakup coincidence cross 
section on 197 Au at Ed = 12 MeV (Fig. 4 of [8*7]). The postacceleration effect can clearly be seen, as 
the maximum of the proton energy (~ 7.5 MeV) is larger than the one of the neutron (~ 2.5 MeV). 
The experimental data are taken from [75] . 

Of course this picture is based on the purely classical interpretation of this process, and will be modified 
in a quantal treatment, where such a breakup point does not exist. The correct semiclassical limit of 
the theory in this case can be found, e.g., in [51] . A purely classical formula for this postacceleration, 
where the breakup point corresponds to the distance of closest approach — i.e., b = r in Fig. |S] — is 
given in j^SJ. Postacceleration is clearly observed in low energy deuteron breakup, in the theoretical 
calculations, as well as, in the corresponding experiments, see Fig. and also, e.g., [8lfl IHTj . 

The formula Eq. (J85|) is also useful for the description of the Coulomb dissociation of halo nuclei at 
high beam energies, see [HE] ■ Within this theory postacceleration effects become negligibly small in the 
high energy region. This is seen in the numerical calculations [88J and in the analytical investigations to 
be described below. It can, e.g., be applied to n Be and 19 C Coulomb dissociation experiments |6lH IBT] 
(We disregard here the importance of finite range effects). 

On the other hand the 1 st order semiclassical Coulomb excitation theory was widely applied in the 
past years to the Coulomb dissociation of high energy neutron halo nuclei, see, e.g., jHH]- The theory 
corresponds to the prior form, mentioned above. The question of higher order electromagnetic effects 
was studied recently in within this framework. These effects were found to be small, for zero range, 
as well as, finite range wave functions of the a = (c + n)-system. Higher order effects were recently 
studied in a post-form DWBA approach [HEUm], which was already discussed in Sec. |3] 




Ep(MeV) 
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We now show that the full CWBA amplitude is directly related to the Born approximation, Eq. (|7Sjl 
or the semiclassical result, see Eq. (J70J). This is true (with modifications to be discussed below) not 
only for r] a <C 1 but also for arbitrary values of i] a and r\ c . For this to be the case, the beam energy 
must be high (compared to the binding energy Eq) and the two fragment must be scattered to forward 
angles. 

As was discussed above, the plane-wave approximation to the post-form DWBA has some similarity 
to the bremsstrahlung process. Replacing the neutron by a photon the diagrams of Fig. [7| are analogous 
to the bremsstrahlung diagrams in lowest order. The electron-photon vertex corresponds to the (c + 
n) —>■ c + n breakup amplitude, the neutral massless photon to the massive neutron. In the case of 
bremsstrahlung it is well known that for arbitrary values oft], even for 77 ^> 1, one obtains the Born 
approximation result, as long as the scattering is into a narrow cone in the forward direction. This 
leads one to suspect that higher order effects are not very large in the case of high energy Coulomb 
dissociation also, where the fragments are emitted into the forward direction. 

The T-matrix can be evaluated analytically in this model due to the well known Nordsieck formula 
[91j . Using this formula one obtains the T-matrix Eq. in terms of a hypergeometric function F as 
well as its derivative F'. The argument of the hypergeometric function F (and F') is given by [77| I88j: 

^ = 2g 2 (g a g c + g a ■ g c ) - 4(g • q a )(q- g c ) ^ 
(q 2 - 2q ■ q a ){q 2 + 2q ■ q c ) 

We observe that the parameter ( is found to be negative and — (0 ^> 1 for perpendicular momentum 
transfers q±_ ^> 2r] a q\\ (nonadiabatic case) and beam energy large compared to the binding energy. 

This was already noticed in the numerical evaluation of the process: the hypergeometric series 
therefore does not converge and an analytic continuation had to be used. Here we use this fact to 
our advantage and make a linear transformation to get the argument of the hypergeometric function 
close to zero. The transformation we are using leads to the argument of the hypergeometric function 
z — 1/(1 — Co) (Eq. (15.3.7) of [H2])- In this respect our approach differs from the one used in the 
bremsstrahlung case, where a transformation giving an argument close to one is used. Using only the 
lowest order term in the hypergeometric series one obtains after some algebra (up to an overall phase) 



T w 4tt D f coul e 



2 s 



Qa ~ (<fn + <fc) 2 m a q 2 c - (<f n - q a y 



Hereby, the relative phase is <fi = (t (t] c ) — o- Q (r] a ) — cr (£) — £/2 log \( \. The 00(77) = arg T(l + irf) are 
the usual Coulomb phase shifts, and £ = r\ c — r\a ■ 

The correspondence to the Born result Eq. (J76|) is clearly seen. One only has an additional prefactor 
e _ 2^ and a relative phase e ±l< ^ between the two terms. The phase is 0(£). Since v c ~ v a the quantity 
£ is usually very small and so is for the cases of, e.g., [HSUE!!- The prefactor is also well known in the 
semiclassical theory, where it accounts for the replacement of the Coulomb bended trajectories with the 
straight line trajectories. Both corrections vanish in the limit £ — > and the result coincides with the 
usual Born approximation (even if rj a and r] c are not small). 

Above we have established the correspondence of the full CWBA result with the (Born) plane wave 
result and with the semiclassical result in the sudden limit. One may expect also to find a connection 
between the two theories in the adiabatic limit, corresponding to q± < 2r] a q\\. 

For i] a , r] c ^> 1 one can define a classical path for both a and c in the initial and final state. (For 
high beam energies, the straight line approximation can be applied and the impact parameter b is given 
by b = 2r] a /q±. The adiabatic case corresponds to the case of large impact parameters b > 1/gii.) For 

lr rhe parameter £ is the one used commonly in NR Coulomb excitation. The parameter £ as introduced above (see 
Eq. 1)11) fl hi Sec. [2| corresponds to the parameter £(#), see, e.g., p. 6 of PP in the literature on NR Coulomb excitation. 
Instead of the scattering angle 9 we use the corresponding impact parameter b. 
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the adiabatic regime (q± < 2rj a q\\) the well known exponential decrease with the adiabaticity parameter 
is observed in the numerical calculations. In this case, the inequality — £ ^> 1 is no longer generally 
satisfied. 

The semiclassical limit can be obtained with analytical methods in this case also. Such a method is 
valid in both the adiabatic and nonadiabatic case as long as r] a ,i] c ^> 1. 

We use two methods to show this result: in the first approach we use the confluence 



together with some relations of to relate hypergeometric functions to (modified) Bessel functions, 
which appears also in the semiclassical limit. 

A second approach makes use of the partial wave expansion of the neutron plane wave. The semi- 
classical limit is made by approximating the radial wave function with the WKB approximation. In the 
case of the plane wave expansion the resulting sum over I can then by done and an analytic expression 
can be found for the absolute value of the matrixelement T [93.. 

A similar situation is encountered in the theory of bremsstrahlung and Coulomb excitation, see 
Sec. II E of j2]. There a fully quantal expression for the differential cross-section for dipole Coulomb 
excitation is given in II E.62. It looks similar to the corresponding expression for Coulomb breakup 
(see jZZI)- The semiclassical version of this formula is found in II E.57. It is noted there that it can be 
obtained from the quantal expression by letting both i] a and r] c go to infinity and perform a confluence 
in the hypergeometric functions. 

Postacceleration effects are also important for Coulomb dissociation studies of radiative capture 
reactions of astrophysical interest. We expect that our present investigations will shed light on questions 
of postacceleration and higher order effects in these cases also. Postacceleration was studied numerically 
in as expected these effects turned out to be small at high energies. 

We also note the following important physics point: in the post-form DWBA we can introduce in a 
straightforward way neutron-target interactions [HI]. One may say that breakup is a kind of "transfer 
into the continuum" , see also On the other hand, in the prior form DWBA the final state 

wave function is the product of a cm. motion and an internal n-core wave function, governed by the 
interaction V nc . In this case the breakup process may be viewed as an inelastic excitation of the bound 
(c + n)-system into the continuum. 

5 Nuclear Effects, Coulomb-Nuclear Interference 

A common problem for the Coulomb dissociation method is the influence of the strong interaction 
between projectile and target. For heavy target nuclei Coulomb excitation often dominates due to the 
coherence factor Z 2 , whereas the total nuclear cross section increases with R 2 ~ A 2 ! 3 , direct reactions 
(elastic nuclear breakup is a special case) scale with R ~ A 1 / 3 , see, e.g., jHZl, as they take place in a 
ring area around the nucleus. 

A good method to avoid the nuclear interaction altogether is to use beam energies below the Coulomb 
barrier, which suppresses nuclear effects very effectively. As was discussed above in Sec.^ one can access 
only low energy nuclear states in this way. We are mainly interested in higher lying states and higher 
beam energies are necessary. Again, we can suppress nuclear effects, this time by going to forward 
scattering angles. This corresponds, in a semiclassical picture, to trajectories where the nuclei do not 
touch each other. The (total) electromagnetic excitation probability falls off at least with 1/b 2 ; it is 
maximal for grazing impact parameters, where also the nuclear effects set in. For too large impact 
parameters, the adiabatic cutoff sets in and the excitation probability drops exponentially to zero. I.e. 
we have the same problem as Wilhelm Tell [US], who had the task to hit a rather well defined spot, with 
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grave consequences in the case of failure. By the wave nature of the projectile nuclear effects cannot be 
totally avoided, they will enter somehow, but often this effect is negligible. Even for small scattering 
angles there can be some effect coming from the interference of that part of the wave function that is 
bent by the Coulomb interaction with that part of the wave function, which is diffracted by the "black 
disc" of the target. This will lead to diffraction patterns, see the discussion in jHHJ. 

It should be kept in mind that one does understand nuclear effects, at least in a semiquantitative 
way. It is certainly best to choose the experimental conditions (beam energy, scattering angle,. . . ) in 
such a way that these effects are virtually negligible. Such considerations can be based on the good 
knowledge of nuclear effects. As a second best approach they can be taken into account with confidence 
by using the results of some kind of nuclear reaction theory. Direct (nuclear) reactions have been 
extensively studied in the past decades and it can safely be said that one understands rather well their 
main characteristic features, but not as well as Coulomb excitation itself. This section is devoted to 
these nuclear processes, Coulomb-nuclear interference and the theoretical methods used to describe 
them. These methods sometimes involve heavy computing. Entire conferences are devoted to this 
subject jlOOj . which we only briefly summarize here. 

5.1 Various Kinds of Nuclear Reaction Mechanisms 

Nuclear breakup is dominant for light target nuclei and is used by itself as a tool to explore the 
structure of exotic nuclei, especially nuclear halo systems. In a nuclear projectile-target interaction 
all of the fragments (typically two) can emerge while the nucleus remains in the ground state. This 
is called diffractive breakup and it is a reaction mechanism which leads to the same final state as 
Coulomb dissociation. In addition to diffractive breakup (diffractive dissociation) there is also stripping 
and absorption. In stripping reactions, only part of the fragments comes out of the reaction in the very 
forward direction, while the rest is absorbed; also all fragments can interact violently with the target 
and are absorbed, see |1011 I1U2| 1103] for details. 

Stripping reactions have been the "workhorse" of the study of halo systems in nuclear reactions, 
mainly due to the fact that they can be interpreted in a simple (sometimes too simple) way: in the 
so-called transparent limit the longitudinal momentum distribution of the core after stripping is related 
to the momentum space probability distribution of the halo nucleon. This picture was more refined in 
recent years [75J E7| lll)4| in order to incorporate the so-called core-shadowing effects. The longitudinal 
momentum distribution is related to the momentum- distribution at the surface,, that is, to the Wigner 
transform at the surface |lU5j . 

Diffractive breakup leads to the same final state as Coulomb excitation. As all fragments are 
generally measured in order to reconstruct, e.g., their relative excitation energy, it is only this process 
which needs to be considered as a background to Coulomb excitation. We concentrate in the following 
on this process. We refer here to a number of reviews 1061 1107| I108| I1U9| IllUj of direct nuclear reactions 
with exotic nuclei. 



Nuclear effects can be treated in complete analogy to the DWBA approach explored for Coulomb 
excitation above. The pure-Coulomb wave function \1/ is replaced by the one where an optical potential 
between projectile and target acts in addition to the Coulomb potential. In addition one has to add 
the nuclear excitation, treated in first order here, for higher order effects see also the discussion below 
in Sec. 15.41 

The T-matrix for inelastic scattering in the DWBA is given by, see, e.g., Eq. (7) of jlllj : 



5.2 DWBA 
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Figure 10: The potential V[R,£) is given by the sum of the (Coulomb and nuclear) interactions between 
the individual components of the target with the projectile Vi(R — 

The cm. distorted waves functions are denoted by The nuclear form-factor Ffi(R) is given by 

F fi (R) = jdt V(R,0 0,(0, (90) 

where £ denotes the internal coordinates, see Fig. EH The initial and final internal wave functions are 
4>ij(0- The projectile-target interaction is given by 

F(iu) = EW*-S- ( 91 ) 

i 

We see that nuclear effects enter in two ways: first there is an influence on the distorted waves of 
the cm. motion due to the nuclear interaction. Second there is an additional breakup component 
due to the nuclear form-factor F. This form-factor can be expanded in terms of multipoles and added 
to the Coulomb multipoles. 

At lower energies the interactions of protons, neutrons, and also light ions with different target 
nuclei have been analysed and tabulated |112j , so that they are usually well known. At higher energies 
(few 100 MeV) the knowledge of the phenomenological optical potentials is rather sparse. Often global 
parametrisations (fitting optical potentials for protons and neutrons) are used for a whole range of 
energies and targets, see, e.g., j!13j . or for individual target nuclei, see |114j for a recent reference. The 
interaction between the projectile nucleus and the target (or some fragments with the target) are then 
calculated within a folding approach of the projectile density with the proton-target and neutron-target 
optical potentials. 

V(r) = J d 3 r lPl (n)V NA (\r-ri\), (92) 

where Vna is the nucleon-target potential. 

Another way is to use a double folding approach over both densities (projectile and target). The 
nuclear interaction is then given by 

V{r) = J d 3 r 1 d 3 r 2Pl {r 1 )p 2 {r 2 )V NN {\f- n + f 2 |), (93) 

where Vnn is the nucleon-nucleon potential. 

Only the nucleon-nucleon interaction Vnn is needed in addition to the nuclear densities p^ 2 - A 
typical example of this approach is the DDM3Y interactions |115t I116j , based on a fit of the nucleon- 
nucleon interaction in nuclear matter |117j together with an assumed density dependence. Only the real 
(elastic) part of the potential can be found in this way and the imaginary part needs to be adjusted. In 
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many cases it is chosen for simplicity to be proportional to the real part, with a proportionality factor 
in the range of 0.5-0.9. This proportionality is in general only a first approximation; especially at lower 
energies there are corrections, which can be described by a polarisation potential, see e.g., |118j . 



5.3 Eikonal DWBA 

At high energies the partial wave decomposition of the wave function of projectile and target is 
rather cumbersome. High energy methods, like the eikonal approach ("Glauber method"), are more 
useful and at the same time accurate. The wave nature of the projectile is still taken fully into account 
and leads to characteristic diffraction effects. The use of eikonal wave functions for the initial and final 
states in the DWBA matrixelement is discussed in (3K|, where also its relation to the semiclassical case 
is explored. Higher order effects can be treated within a coupled channels method jSlj, a corresponding 
computer program is described there. In this approach the differential cross section can be written as 
(see also Eq. (|27jl): 

ftnei = «'* / bdbJ^qb) exp(i X (&))a„(&). (94) 

J 

The index // denotes the angular momentum transfer to the target in the beam direction. The excitation 
amplitudes a^(b) are calculated in the semiclassical straight line approximation for impact parameter 
b. The eikonal phase x(b) = Xcip) + Xn{P) takes Coulomb, as well as, nuclear effects into account. It 
involves the Coulomb phase given by 



Xc (b) = -2r]K (kb) « 2r)ln{kb) (95) 
and a nuclear phase, which can be calculated from an optical model potential U opt {r) as 

Xn(P) = -i- r dzU opt (r = VWT^). (96) 
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The simplest way to take nuclear effects into account is the "black disk model", where one sets 
expix = for b < R\ + R 2 and x — Xc for b > R± + R 2 . This is justified due to the strong imaginary 
part of the nuclear interaction at these energies. 

Diffraction effects in this black disk model are studied numerically in |3S] and |37] . In Figs. 1 and 2 
of jHZj the importance of nuclear diffraction can be seen. A reduced integral is defined there in Eq. (8); 
it depends on the characteristic parameters 77 and £. Plots are given for a reduced scattering angle 9/9 gr 
where 9 gr = 2rj9diff and the diffraction angle is given by 9 = l/(kR). It can be seen in theses figures 
that for r) 3> 1 the semiclassical approximation is very good, especially for small values of £. 

A model which takes into account the smoothness of the surface is discussed in (22], where surface 
normalized Gaussians are used for the density of the nuclei. In this way analytic expressions for the 
eikonal phase can be found of the form 

X (b) = xo exp [-b 2 /(al + a 2 2 )] , (97) 

with 

7i 2 v NN p 1 (0)p 2 {0)a 3 1 a 3 2 

Xo = in , 2 , 2 ^ 98 

10(af + a 2 ) 

and di and p«(0) are adjusted to reproduce the experimental nuclear density at the surface of the nucleus; 
for a tabulation see (22] • cttvat denotes the average nucleon-nucleon cross section. 

An advantage of the (Glauber) eikonal approach is the fact, that one does not need a nuclear 
potential in order to derive the eikonal phase. Scattering at high energies is concentrated in the forward 



31 



direction. A knowledge of the elastic scattering amplitude f e i{q) is sufficient, as it is related via a Fourier 
transform to the eikonal phase 



if -* 
exp [ix(b)\ - 1 = y~l J d2 Qfei(Q) exp(-igb). (99) 

In order to determine the eikonal phase in nucleus-nucleus scattering xaa the so called tpp formalism 
can be used, where the interaction is described through a folding of the densities together with the 
nucleon-nucleon scattering amplitude |119j . cf. to Eq. (|93[). In terms of the nucleon-nucleon profile 
function Y NN {b) = 1 — S^N^b) = 1 — exp [iXNN{b)} it can be written as 

XAA(b) = i J d i r 1 d z r 2 p{r l )p{r 2 )Y NN (r^ l - f ±2 - b), (100) 

This is related to the double-folding approach, see Eq. (}9*3*|) . In the limiting case of a zero-range 
interaction one gets 

T NN (x) = 5 (2) (f)-(l - ia NN )a NN , (101) 

where o~nn is the total nuclear cross section and onn the ratio of the real to imaginary part of the 
forward scattering amplitude. Parameterizations of this forward amplitude can be found in the literature 

The tpp formalism, which is often also called the optical limit approximation can be derived from 
the multiple particle Glauber formalism by taking the first term of the cumulant expansion. One wants 
to find an eikonal Xaa, which reproduces the elastic scattering amplitude of the microscopic multiple 
scattering model, see Eqs. (jl03j) and (jl04j) below, 



exp 



iXAA(b)} =<0|n$(&i)|0>. (102) 



By taking the logarithm of both sides and expanding to first order in powers of T^n, see also |121j . 
one obtains Eq. (jlOOJI above. In this way the interaction is approximately included in all orders. A 
discussion of this approach can be found also in \V2'2\ 11231 H05j . 



5.4 Higher Order Effects 

Whereas the Coulomb interaction is long ranged and often the major contribution comes from larger 
impact parameters, where the interaction is not strong and first order theories are adequate, the nuclear 
interaction is short ranged and strong. Therefore higher order effects can become important and should 
be taken into account. 

Within the eikonal DWBA higher order effects in F(R) can be taken into account as discussed above, 
e.g., in the (semiclassical) coupled channel approach as given in j33]. In the case of the deuteron, 6 Li and 
other halo systems, quite elaborate calculations have also been made in terms of the CDCC (continuum 
discretized coupled channel) see |124| I125j , or [1261 11271 H28j . In both approaches one needs to identify 
the dominant structures especially in the continuum in order to restrict oneself to a manageable number 
of channels to be incorporated in the calculation. 

A theory which is well suited for intermediate and high energy nuclear breakup of halo nuclei is 
the Glauber multiple scattering theory [12111122] . In this approach one uses the eikonal approximation 
together with the sudden limit to describe the elastic, as well as, inelastic scattering from the ground 
state to a state m as 

fmo{q) = ^~J d 2 bex P (iqb) J d^miOMO (S(b,S ± ) - 1) , (103) 
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Figure 11: The "wounded nucleus" picture of nuclear elastic breakup (diffractive breakup) of a halo 
nucleus. Due to the strong absorptive (imaginary) part of the nuclear interaction, those parts of the 
wave function, which are touched by the target nucleus are set to zero. This means that the nucleus 
is no longer in its ground state but also in excited states. The projection of this excited state to the 
continuum states gives the amplitude for diffractive dissociation, see also Eq. (jl(J5|) . 



where k is the momentum of the projectile and q = ki — kf. The profile function S depends on the 
individual impact parameters hi = b — It is given by the product of the individual profile functions 

S(b,Z ± ) =pft) = IIexp(%(6i)). (104) 

i i 

It is possible to evaluate this expression numerically with fully microscopic wave function <f) and <ft m , 
see e.g. [T23] . 

An intermediate model, which is well suited for halo nuclei, is the Serber model |lUlj . improved 
by Glauber |lU2j . see also |lU3j : in this model only the relevant clusters are taken into account in the 
Glauber multiple scattering model. 

In the eikonal DWBA, see Eq. a global projectile-target profile function is used to describe the 
cm. motion of the projectile. In the improved Serber model, see Eq. (|1U5|) below, one uses individual 
profile functions S^&j) for each fragment i, e.g., the core c and the halo neutron n \YdO\\l'dl\\l'6'2\\YS'S\ W7\. 
The individual profile functions can then be calculated in the same way as in Eqs. (JHSj), (J57)l or (jlOOJ) . 
The differential cross section for the elastic breakup {diffraction) is then found to be, see also Eq. ()103J) . 



ffiil) = 7T- d bexp(iqb) (f\ exp (iXc(b c )) exp (%(&„)) - 1 |i) , (105) 
2m J 

where b c = b + m n / (m n + m c )r± and b n = b — m c j (m n + m c )r± are the individual impact parameter 
of the core and the neutron with the target, respectively. In this way the most important parts of the 
internal structure of the projectile nucleus are taken into account. This approach can easily be extended 
also to the case of more than two clusters |134t 1135] . Only initial and final state wave functions but no 
intermediate states are needed in this way. 

A simple picture of the diffractive breakup can be found in the limit of the "black disc" model, 
where the individual profile functions are approximated by black discs, see Fig. [TTJ 



5.5 Coulomb- Nuclear Interference 

Nuclear effects enter most sensitively through the interference of the nuclear and Coulomb amplitudes. 
This has been addressed in the past in a number of papers, some recent references are |136H137tll38l fl39j. 

A simple (analytic) result is given in [7j for a neutron halo nucleus. The amplitude of Coulomb 
excitation in first order (Eq. (3.4.16) of [7j) is combined with the nuclear breakup amplitude of Akhiezer 
and Sitenko |14Uj (Eq. (3.4.14) of [7j). In this way a good estimate of the interference effect can be 
found. 
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Figure 12: Fingerprints of Coulomb- nuclear interference: effects of nuclear and Coulomb-interference are 
shown for the breakup of 11 Be on 208 Pb at 72 AMeV. Shown are the contribution from pure Coulomb 
(dashed), nuclear (dotted) and their incoherent (dashed-dotted) and coherent sum (solid line). The 
magnitude of negative and positive Coulomb-nuclear interference terms are shown as inverted triangles 
and plus signs, respectively. The experimental results are taken from jHZj. This figure is adapted from 
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Figure 13: Effects of nuclear and Coulomb-nuclear interference are shown for the breakup of Be on 
208 Pb at 72 AMeV in the projectile (left)and target (right) frame. Shown are the results for nuclear 
breakup (dot-dashed), Coulomb (long-dashed) and Coulomb-nuclear interference (dotted) and the sum 
(solid line). The experimental results are taken from jHZ]. This figure is taken from |13bJ . Fig. 4 there. 
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More detailed calculations have been done in the mean time. E.g., in |136| 1141] a calculation within 
Glauber theory is presented. In [138 the time-dependent Schrodinger equation in the semiclassical 
approximation is solved numerically and compared to a first order calculation and to the eikonal ap- 
proximation. 

In the case of halo nuclei the condition of no nuclear penetration may not be well fulfilled and the 
decomposition of the Coulomb interaction into r< and r> parts (see Eq. (j3J)) may not be reliable, as 
shown in |138j for the case of 8 B. 

We illustrate the fingerprints of Coulomb-nuclear interference and nuclear effects from two recent 
theoretical model calculations. The results shown in Fig. ^] are based on a post-form DWBA approach, 
see Sec. |U The results shown in Fig. [Tulare based on a Glauber type calculation, see Eq. (JTU3J), taking 
higher orders in the nuclear interaction together with the Coulomb interaction in lowest order into 
account. 

In both cases breakup of 11 Be on a 208 Pb target at 72 AMeV are shown in comparison to the 
experimental data of jUZj- One sees that both nuclear and Coulomb-nuclear effects are rather small in 
this case. 

6 Coulomb Dissociation and Nuclear Structure 

6.1 Primakoff Effect 

We want to recall first at this point that the strong nuclear Coulomb field has also been used in particle 
physics to study photon interactions with (unstable) particles. Typically a high energy secondary beam, 
like a beam of A particles hits a heavy target nucleus. In the nuclear Coulomb field, high energy E° 
hyperons are produced in the reaction 

A + Z -> S° + Z. (106) 

The cross section for this process is proportional to the B (Ml)- value for the A —>■ S° electromagnetic 
excitation. From the measurement of this cross section the lifetime of the S° is obtained. For further 
details of the analysis see [7j, where further references and other applications of the Primakoff effect (like 
the determination of the pion polarizability) are given. A more recent example is the radiative decay 
width measurement of neutral Kaon excitations using the Primakoff effect by the KTeV Collaboration 
[142 . They used Kl mesons in the 100-200 GeV energy range to produce the axial vector (1 + ) mesons 
_Ki(1270) and 1^(1400) in the nuclear Coulomb field of a Pb target. In this way the radiative widths for 
the decay of these particles into K° + 7 could be determined. For further details we refer to this paper 
|142j . This approach is very similar to the one used nowadays in nuclear physics: an exotic (unstable) 
nucleus is excited by the quasireal photons provided by the nuclear Coulomb field. 

We mention a second example to illustrate the use of equivalent photons: At Fermilab the Primakoff 
effect is used for a determination of the proton polarization of a 185 GeV/c proton beam by means of 
azimuthal asymmetry measurements |143j . This method could also be useful for antiproton beams. For 
further details we refer to this reference (143] . From this work one can see that it is possible to study the 
photon-proton interaction in the nucleon resonance region using equivalent photons. In Fig. El w e show 
the invariant mass spectrum of the tt° — p system in the reaction p + Pb — > n° + p + Pb for a momentum 
transfer |i| < 1 x 10~ 3 (GeV/c) 2 . One can clearly recognize the nucleon-resonances (especially the A- 
resonance) which are excited with the continuous equivalent photon spectrum provided by the Coulomb 
field of the Pb target nucleus. Since the Lorentz factor 7 of the proton is about 200 in this experiment, 
the corresponding equivalent photon spectrum is quite hard: (see the discussion in Sec. 12.21 about the 
cutoff in the photon spectrum at u> max ~ j/R) we can roughly take 1/R ~ 30 MeV and obtain a 
maximum photon energy of about 6 GeV. 
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Figure 14: The invariant mass spectrum of the 7r° — p system in the reaction p + Pb — > 7r° + p + Pb 
with a 185 GeV/c proton for a momentum transfer of \t\ < 1 x 10 _3 (GeV/c) 2 is shown. These small 
momentum transfer values ensure that photon exchange is the dominant excitation mechanism. Peaks 
due to the A + (1232) and iV*(1520) resonances are clearly seen. Reproduced from Fig. 2 of [143], to 
which we refer for further details. Copyright (1990) by the American Physical Society. 

6.2 Some Aspects of Electromagnetic Excitation in Relativistic Heavy Ion Collisions 

Electromagnetic excitation is also used at (ultra-) relativistic heavy ion accelerators to obtain nuclear 
structure information. Recent examples are the Coulomb fission studies of radioactive nuclei at GSI 
|144[ 1145] and Coulomb fission of 208 Pb j!46j at SPS/CERN. Due to the logarithmic rise of the cross- 
section with beam energy, cross-sections for the excitation of the giant dipole resonance ( Weizsacker- 
Williams process) at the relativistic heavy ion colliders RHIC and the forthcoming LHC(Pb-Pb) at 
CERN are huge [HH EHl EH] , of the order of 100 b for heavy systems (Au-Au or Pb-Pb). 

The neutrons from GDR decay were observed at RHIC (150] . In this reference mutual Coulomb 
dissociation was measured in ^snn = 130 GeV Au-Au collisions. This is shown in Fig. (taken from 

urn)- 

The two classes of events are defined according to the number ubbc of hits in the BBC [beam beam 
counters) photomultipliers. The PHENIX (one of the detectors at RHIC) BBC measures the relativistic 
charged particles produced in cones around each beam with a rapidity range 3.05 < \r}\ < 3.85) with 2n 
azimuthal coverage. For hadronic events we have Ubbc > 1 in each arm, for Coulomb events there is 
t^bbc < 1 in at least one arm. We see that the "Coulomb" events tend to have a low neutron multiplicity. 
This is easily understood qualitatively: Coulomb events are mainly due to the electromagnetic excitation 
of the GDR (in 197 Au) in ultraperipheral collisions (UPC). The GDR decays subsequently by the 
emission of one (or only a few) neutrons. See |150j for further details. 

In colliders, this effect leads to a loss of the beam, due to the particle (neutron) decay of the GDR. 
Even worse, this effect can lead to a localized beam pipe heating in Pb-Pb collisions at LHC, as was 
noticed in |151j . (An even more severe process in this context, limiting the maximum luminosity that 
can be achieved at the LHC, is the electromagnetic process of bound-free electron pair production.) On 
the other hand, this effect can also be useful as a luminosity monitor by detecting the neutrons in the 
forward direction as first proposed in |152j and demonstrated in (150] . One measures the neutrons which 
are produced in the decay of the giant dipole resonance, which is excited in each of the ions (mutual 
excitation). Since this process has a steeper impact parameter dependence than the single excitation 
cross-section, there is more sensitivity to the cut-off radius and to nuclear effects. The neutrons from 
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Figure 15: Single arm ZDC (Zero Degree Calorimeter) neutron multiplicity spectrum for Coulomb and 
hadronic events. The two types of neutron emission can be separated from each other as hadronic 
production also leads to charged particle production in the more central parts of the detector. See text 
and |15Uj for details. Reproduced from Fig. 2 of |15Uj . Copyright (2001) by the American Physical 
Society. 

the GDR decay can also serve as a trigger on the ultraperipheral collisions. For details and further 
references, see |1481 1149] . 

UPCs are also of great interest in particle and hadron physics: at RHIC and LHC the equivalent 
photon spectrum extends to about 500 GeV and 1 PeV (10 3 TeV), respectively, in the rest frame of one 
of the ions. This leads to many interesting applications like the coherent production of vector mesons 
|153[ 1154] . Coherent p° production with and without simultaneous GDR excitation has been observed 
at RHIC |155j . This is a very hot topic, but it is outside the scope of the present article, we refer to 
the reviews in [HH EH EH [T37| . 

6.3 Intermediate Energy Coulomb Excitation of Discrete Levels and Gamma Decay 

Electromagnetic excitation is a very clean and efficient way to excite nuclear states. Interesting nuclear 
structure properties like the position of energy levels, as well as, electromagnetic transition matrix 
elements can be determined. The general features of electromagnetic excitation have been given in 
Sec. 121 Collective states, i.e., those with large electromagnetic matrixelements, are most strongly 
excited. While in Coulomb excitation below the barrier one can only populate low lying states like 
rotational or low lying vibrational states (see discussion in Sec. 12}, it becomes possible at intermediate 
energies to excite also high lying states, the prime example being the GDR. The main features are 
shown in Fig. El taken from see also Figs. 3.3 and 3.7 of jjj- In these figures the excitation cross 
section of an exotic nucleus with typical values for structure parameters (whose exact values do not 
matter now) is shown. The low lying 2 + -collective level is excited most strongly at the low energies, 
while the GDR takes over at the higher energies. Two- and three-phonon GDR's (labelled 2-ph GDR, 
and 3-ph GDR) are excited with smaller cross-sections (see Sec. I3.ip . They tend to a constant in the 
limit of high beam energies. The one-photon GDR excitation cross section first rises quite steeply with 
increasing energy, overcoming the adiabatic cut-off criterion; eventually it will show the asymptotic rise 
with In 7, see Eq. ()21|). The (high lying) isoscalar quadrupole excitation is also given. The cross section 
for the low lying 2 + excitation shows a l/w 2 -behaviour, while its logarithmic rise will only occur at 
extremely high energies (see the remarks in Sec. EI). These excited states can decay by photon- (or, at 
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Figure 16: Cross sections for the electromagnetic excitation of collective states of an exotic nucleus on a 
Pb target. Typical nuclear structure parameters for a medium mass nucleus were adopted. This figure 
is taken from jSj, p. 128. 

sufficiently high excitation energies, by particle-) emission. 

Photons emitted in-flight from intermediate energy projectiles can be distinguished with modern 
gamma ray detectors from 7-rays originating from the target. In order to accomplish this, the in- 
formation on the emission angle is used to correct for the Doppler broadening. The development of 
position-sensitive 7-detectors was a key to the success of this method. For a detailed description of this 
field we refer to the review by Glasmacher |158j and more recently [159 and j^j. 

This method is very well suited for radioactive beams and very valuable information about the 
nuclear structure of unstable nuclei has been obtained and is expected to be obtained in the future. 

Many interesting examples of the investigation of nuclear structure effects of unstable nuclei with 
this method are given in [159 . We refer to this reference for more details. Let us mention here only 
some points: 

Electromagnetic excitation of intermediate energy (exotic) beams has been developed into a useful 
spectroscopic tool [16QJ I161j . This method is ideal to study the behaviour and onset of collectivity 
(deformation) for + — > 2 + transitions in nuclei far from stability. By measuring the excitation energies 
of the first 2 + states and the corresponding B(E2)-va\ues, nuclear structure effects like deformation, 
can be studied in a unique way for nuclei far off stability. 

Electromagnetic excitation of the I s * excited state in 11 Be has been studied experimentally at GANIL 
[162 , RIKEN |163j and MSU |164j . This is a good test case, since the B(El)-vahie of the corresponding 
ground-state transition is known independently. Theoretical calculations [165] I166j show that higher 
order effects are small. They decrease with increasing beam energy, as expected. 

Subsequently, many interesting nuclear structure studies followed at RIKEN, GANIL, MSU/NSCL 
and GSI with this method, which can be considered as well established by now. 

We mention here that the neutron-rich nucleus 34 Mg has been studied recently via Coulomb excita- 
tion using a radioactive beam on a Pb target |167j . The high B(E2) value for the excitation of the 2f 
state corresponds to a quadrupole deformation parameter fa of 0.58(6), implying an anomalously large 
deformation of 34 Mg. 
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The work done at MSU can be found at |168j . Let us mention one of their latest results: The 
0+ s — > 2 + excitations in the mirror nuclei 32 Ar and 32 Si were compared to each other. This is a sensitive 
test of isospin symmetry which could be extended to a T = 2 isobaric multiplet |169j . 

In this context it is worth mentioning that a computer program for scattering at intermediate energies 
has recently been published in [31]. Nuclear, as well as, Coulomb excitation processes are calculated 
using the eikonal approximation and semiclassical coupled channels methods. The eikonal method is 
very appropriate for intermediate energies: it is both numerically accurate and fast (as opposed to 
DWBA methods, where also the relative center of mass motion, with its huge number of partial waves, 
is treated in a quantum- mechanical way). The general problem of the coupled channels method is to 
find the important states to be included. In many practical cases this proves to be possible. In any 
case, as we have shown above, the strength parameter is proportional to 1/v, i.e., higher order effects 
tend to be small for intermediate and high beam energies. On the other hand this also means that 
one cannot have efficient multiple Coulomb excitation for intermediate beam energies. The two-photon 
excitation of the strongly collective DGDR discussed above is an interesting exception. For low energy 
Coulomb excitation, the multiple excitation of rotational bands is a very well known feature. Angular 
distributions of elastically and inelastically scattered particles, as well as, angular distributions of 7-rays 
are also calculated in [3"4] . 



6.4 Intermediate Energy Coulomb Dissociation, Invariant Mass Spectroscopy and 
Low Lying El Strength 

Nuclei in the valley of stability show a prominent collective mode, the Giant Dipole Resonance (GDR). 
In this mode, all protons oscillate against all neutrons. This state has an excitation energy of about 
Egdr ~ 80MeV A' 1 / 3 and exhausts the classical energy-weighted Thomas-Reiche-Kuhn (TRK) sum 
rule to a large extent. As a consequence of this, low lying .El-strength is strongly hindered. In |170j it 
was shown that this hindrance of El-transition strength disappears when loosely bound nucleons are 
involved. This is, e.g., the case in neutron halo nuclei. A classic case is the El-strength in the prototype 
of a neutron halo nucleus: the deuteron. Using simple zero range wave functions, the El strength is 
calculated in |171j . A similar calculation for neutron halo nuclei is given in [7j, which clearly shows how 
the low-lying El-strength is obtained in a single particle model. With the zero range bound state wave 
function, see Eq. (JSHJ) above or Eq. (3.4.11) of [7j 




(107) 



and the continuum wave function, see Eq. (}62|) or Eq. (3.4.12) of [Zj, (for the relevant partial waves 
(I > 0) this corresponds to the spherical Bessel functions, which describe a free particle) we obtain for 
the electromagnetic matrix-element Eq. (see Eq. (3.4.19) of [7 ): 

M(Elm) = eV2^{-i) l l\2 l+ \Z^[ + (-1)%/^) . g ' Y lm (q), (108) 

(k z + q A ) t+i 

with (3i = m%j (mi + m?) an d $2 = mi/ (mi +7712) and Z\ and Z 2 denote the charge numbers of the two 
clusters. 

A general discussion of dipole strength in neutron rich nuclei is given by Hansen and Jonson [172 
and Ikeda |173j . A significant low- lying El strength is found. These modes are sometimes called soft 
dipole state or "pigmy resonance", see Fig. El An early discussion of the disappearance of the hindrance 
of El-transitions involving loosely bound nucleons is given in |170j . 

A measure of the El-strength of certain configurations is given by the cluster sum rule: in the usual 
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Figure 17: In the GDR all neutrons move against all protons (hard mode, lower figure). In a nucleus 
which contains also loosely bound nucleons, e.g., in a 2-neutron halo nucleus like n Li, a new type of 
oscillation, called the soft dipole mode (upper figure), can occur. In this mode the two loosely bound 
neutrons oscillate against the inert core of strongly bound protons and neutrons. This gives rise to El 
strength at low energies. Figure adapted from [173 . 

TRK sum-rule the relevant degrees of freedom of the nucleus are the neutrons and protons 2 . In this 
case the TRK sum rule tells us that 

S^-^i'-Stt- (109 » 

where N ,Z, and A denote the neutron, proton and mass number of the nucleus and where D n0 is the 
dipole matrix-element between the ground state and the excited state n. In this equation we have 
neglected the effect of exchange forces in the nucleon-nucleon interaction, which is of the order of some 
10%. We can determine the contribution of certain cluster configurations to the energy-weighted sum 
rule by a so-called cluster sum rule: assuming an inert core with mass M and charge Z and some 
valence neutrons with mass m as the relevant degrees of freedom, one can write down a cluster sum 
rule, taking those as the relevant degrees of freedom. In this sum rule the factor NZ/A in Eq. ()109|) is 
replaced by ^a^+a^mm ' w h ere an d M = (iVj + Zj) are the neutron, proton, and nucleon number 

of each cluster |175j . In the case of a neutron halo [Zi = 0, either single neutron or more) this gives a 

factor (a^Im)a 1 ■ 

Effective field theories (EFT) are nowadays also used for the description of halo nuclei, sec [176J. 
The relative momentum k of the neutron and the core is indeed much smaller than the inverse range 
of their interaction 1/R and kR is a suitable expansion parameter. (In our model of the pure-Coulomb 
breakup of a bound state bound by a zero range force, see Sec. |U above, we have R = 0, i.e., kR = 
and we have the zero order contribution of the expansion). Effective range theory seems a natural 
starting point. This aspect was pursued in \177\ 1178] . In |178j radiative capture cross sections into s-, 
p- and d-bound states are calculated in simple models, and the cross sections depend only on a few low 
energy parameters. The neutron halo effect on direct neutron capture and photodisintegration of 13 C 
was studied in |179j and [180J. In their figures it can very well be seen that the radial integrals are 
dominated by the outside region. While they find a sensitivity on neutron optical model parameters 
for s — > p capture, this sensitivity is strongly reduced for the p — > s and p — > d capture cases. In 1 1 7(ij 
it is remarked that the EFT approach "is not unrelated to traditional single-particle models" and that 
"it remains to be seen whether these developments will prove to be a significant improvement over more 
traditional approaches. " With a wealth of data on halo nuclei to be expected from the future rare ion 
beams we can be confident that these questions will be answered. An interesting effect is known from 
the deuteron, the mother (prototype) of all halo nuclei: for very low energies, there is an Ml s — ► s 

2 On a more microscopic scale the relevant degrees of freedom are the quarks, see also Exercise 3 of |174) . The nucleon 
resonances, see Fig. ^] are a manifestation of these degrees of freedom 
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Figure 18: Chart of the (particle stable) light nuclei. Stable nuclei, one- and two-neutron halo nuclei, 
four-neutron cluster nuclei and proton halo nuclei are marked. The deuteron (marked in black) is the 
prototype of a one-neutron halo nucleus. Coulomb dissociation experiments of most of those nuclei are 
discussed in the present and the following sections. 

transition, which dominates over El |171j . It remains to be seen whether a similar situation will be 
found in other halo nuclei. 

In Fig. EH we give a part of the chart of nuclides. We concentrate on the light nuclides and indicate 
one- and two- neutron halo nuclei. On the proton rich side, 8 B (and probably also 17 F, at least its l/2 + 
first excited state) can be considered as a proton halo nucleus. 8 He, as well as, 14 Be can be regarded as 
four- neutron cluster nuclei. 

Coulomb dissociation of exotic nuclei is just the right tool to study this nuclear structure problem. 
By this method one determines the electromagnetic matrix elements between the ground state and 
the nuclear continuum. A classic example is the 11 Be Coulomb dissociation. The excitation energy 
spectrum of the 10 Be+ri system in the Coulomb dissociation of the one-neutron halo nucleus n Be on a 
Pb target at 72 AMeV was measured jHZ] , from which we show Fig. la-c as Fig. 

Low lying El-strength was found. It shows the shape of the B(E1) distribution obtained from the 
simple zero range model of Eq. (|1(J8|) . This distribution of low lying El-strength is one of the most 
convincing demonstrations of the neutron halo phenomenon. 

The Coulomb dissociation of the extremely neutron-rich nucleus 19 C was recently studied in a similar 
way [66J. From the shape of the dipole distribution, the neutron separation energy of 19 C could also 
be determined to be 530 ±130 keV. As can be seen from Eq. (|108j) . this shape depends directly on the 
binding energy parameter k. 

We note the similarity of the shape of the low lying B(E1) strength in the case of n Ee (Fig. [TT?j) . 
19 C (Fig. [S] above) and also for the deuteron (see Fig. 4.1, p. 609 in j!71j ). The reason is that all these 
nuclei are halo nuclei: the ground state wave function is an s-wave function with a small binding energy 
parameter k and it is rather well described by Eq. (|107|) or (|58|). The neutron spends most of its time 
in the classically forbidden region. The continuum is a p-wave function which differs little from the 
spherical Bessel function which describes a free neutron. At the low relevant neutron energies there is 
only little interaction with the strong short range force in the p- or higher partial wave channels. In 
this case, the electric transition dipole moment is given by Eq. ()108|) above. 

The El-strength in other neutron-rich carbon isotopes has also been studied experimentally: Cou- 
lomb dissociation of 15 C is again well described by a model with a 14 C-core coupled to an si/ 2 neutron. 
The analysis of |181j gave a spectroscopic factor of 0.75, consistent with the one found from (d,p)- 
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Impact parameter b [fm] 

Figure 19: Figl(a) Dissociation cross sections of n Be on a Pb target as a function of the excitation 
energy E x for the 10 Be+n system. Data from a C target are also indicated, (b) Dipole strengths deduced 
from the daco/dE x spectrum are shown by the open circles, (c) Impact parameter dependence of the 
Coulomb dissociation cross section. The low lying El strength in n Be is clearly seen. This figure is 
taken from [HZj, Fig. 1 there, where further details can be found. 
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Figure 20: The excitation cross section da/dE* is shown for the breakup of 6 He at 240 AMeV on 
a Pb target, where E* has been reconstructed from the invariant mass of the a + n + n fragments. 
The dotted curve shows the result of a semiclassical perturbative calculation using the dB(El)/dE* 
distribution from the three body model of |186| , the solid one shows the convolution of the dotted curve 
with the instrumental response. The two arrow show two 2 + resonances, a known one at an energy of 
E* =1.80 MeV and a predicted one at E* =4.3 MeV [Wj. Reproduced from Fig. 3 of [HE], where 
further details can be found. Copyright (1999) by the American Physical Society. 



reactions. It would be of interest to compare their results obtained with a 606 AMeV 15 C beam to those 
of |182j . see also Sec. 

The Coulomb dissociation of 17 C was also studied. One can adopt a model of a neutron coupled to a 
16 C core. The importance of the 2 + -core excited state is directly evident from the Coulomb dissociation 
measurement: Coulomb dissociation is dominantly accompanied by the emission of photons from the 
1.77 MeV 7-line from the deexcitation of the 2 + state in the 16 C core. A sizeable s-component of the 
neutron coupled to the 2 + core is found, and it is concluded that the ground state spin of 17 C is either 
3/2 + or 5/2 + . For further details we refer to |181j . 

It is appropriate to mention at this point also a different exotic nucleus: the hypernucleus |H. The 
binding energy of the hyperon to the deuteron core is 0.06 ±0.06 MeV. The Coulomb disintegration 
of this hypernucleus depends sensitively on its binding energy, see |183| I184| I185j . 

Quite similarly, the Coulomb dissociation of the 2n-halo nuclei can be studied experimentally in 
order to investigate the low lying .El-strength. Certainly, this is a richer field than the study of 
single-neutron halo nuclei, where essentially single particle effects show up. Now there are in addition 
interesting correlation effects. E.g., 11 Li and also other 2n-nuclei like 6 He are so-called Borromean 
systems: the core-2n-system is bound, whereas none of the binary subsystems (n-n, or core-n) are 
bound. 

The most prominent example is n Li (for 6 He see below), which was studied experimentally in various 
laboratories |191|. 11921 [193J. The momentum distributions for the ( 11 Li, 9 Li + n + n) breakup reaction 
were studied theoretically in ^94 and |195j . In an experiment at MSU |196j . the correlations of the 
outgoing neutrons were studied. Within the limits of experimental accuracy, no correlations were found. 
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Figure 21: Photo-neutron cross sections measured |189j for the unstable oxygen isotopes 20 > 22 O in 
comparison to that of the stable isotope 16 0. The data for 20 > 22 O are compared to shell model calculations 
190 . The (7,p) thresholds are indicated by arrows. This figure is taken from [2j p. 86. 
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(To be sure there must be correlations, as the Borromean effect itself shows, they are just hard to see 
in the breakup experiments). Further studies were reported in 1193] . In this reference invariant mass 
spectroscopy of 10 Li and n Li was performed. It was concluded that "a strong di-neutron correlation can 
be disregarded" . Low lying El- strength was found, it exhausts up to about 8% of the TRK-sum rule 
below 4 MeV excitation energy which corresponds to 96% of the cluster sum rule assuming a 2n-halo. 
Two low lying structures were observed for 10 Li (n- 9 Li final state interaction), see also [197J. 

The Coulomb breakup of the neutron-rich He-isotopes 6 He and 8 He has also been investigated 
experimentally \19S\ I188j . see also Fig. |2"U1 

The Coulomb dissociation cross-section at 227 AMeV on Pb is about a factor of three smaller in 8 He 
than in 6 He. This indicates that 6 He is a halo-nucleus, with low-lying El-strength, whereas 8 He is not. 
The low lying strength in 6 He (up to an excitation energy of 5 MeV) was found to exhaust 10 ± 2% of 
the TRK sum rule, which corresponds to about 40% of the cluster sum rule. 

The dissociation reaction 6 He — > a + In was studied at 23.9 AMeV on C, Al, Cu, Sn and Pb 
targets 199J . For U the Coulomb part accounts for 2/3 of the total two- neutron removal cross section 
°"-2n- Only a small correlation between the two neutrons was found. The 5(i?l)-strength found is in 
agreement with |188j . an experiment with about ten times the beam energy of |199j . 

Similar studies were also carried out in GANIL [200 . The neutron-neutron correlation there was 
analysed in terms of the radius of the "producing source and using Dalitz plots, see also |181j . Final state 
effects are taken into account. Differences between the Coulomb and nuclear breakup are found, which 
are interpreted as being due to the direct breakup in nuclear dissociation and a resonant breakup in 
Coulomb dissociation. 

Photoneutron cross sections for unstable neutron-rich oxygen isotopes were recently studied in 189J. 
Their results for the isotopes 20 > 22 O are shown in Fig. where also the corresponding results for the 
stable isotope 16 are given. It was found that there is systematically a considerable fraction of low-lying 
dipole strength. Low lying El strength was also found in the proton-rich isotope 13 by a 13 — > p+ 12 N 
Coulomb dissociation measurement |201j . It is very promising to extend such types of studies to even 
heavier unstable neutron- (and also proton-) rich nuclei. This is also a challenge for modern nuclear 
structure theory. While we just discussed the single particle aspect, which is the main point for 1-n, or 
2-n halo nuclei, the onset of a pigmy collectivity is a question for more involved theories. The exploration 
of the structure of unstable nuclei is a large field for modern nuclear structure theory, which benefits 
also from the modern high computing power. This field is rapidly expanding, as can be seen from the 
many conferences which are taking place on this topic; as recent ones let us mention the "Hirschegg 
workshop 2003 on Nuclear Structure and Dynamics at the Limits" |202j or "INT workshop on Reaction 
Theory for Nuclei Far From Stability"[100, r We can refer here only to further references which may, 
e.g., be found in j^j. 



7 Nuclear Astrophysics 

Du verstehst nicht die Sterne 

ohne die Kerne. 

(unknown poet, last century) 



In nuclear astrophysics, radiative capture reactions of the type 

b + c-^a + j (110) 
play a very important role. They can also be studied by the time-reversed reaction 

7 + + (111) 
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at least in those cases where the nucleus a is in its ground state. 

Such a photodissociation experiment has been recently performed with a real photon beam by 
Utsunomiya et al. |2()3j . They investigated the photodisintegration of 9 Be with laser-induced Compton 
backscattered 7 rays. We refer to this work for further details. 

The two cross sections Eqs. (|110|) and (Jill}) are related to each other by detailed balance 
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g(t + ^ B + 7) = -^g_ .^ (7 + a ^ t + e) , (112 ) 

where the wave number in the (b + c)-channel is given by k 2 CM = 2fjL bc E C M with fi bc the reduced 
mass; the photon wave number is given by fc 7 = (Ecm + Q)/(hc), where Q is the Q- value of the capture 
reaction Eq. ()11U|) . One typically has k^/kcu *C 1 for energies not immediately in the threshold region 
kcu ~ . Due to this phase space factor there is a strong enhancement of the dissociation cross section 
Eq. (|111|) as compared to the capture cross section Eq. JTTUJ). 

As a photon beam, we use now the equivalent photon spectrum, which is provided by the Coulomb 
field of the target nucleus in the fast peripheral collision There are previous reviews, where both 
experimental, as well as, theoretical aspects have been given, see |T7J HHj. We want to concentrate 
here on theoretical aspects, but also discuss experimental issues, at least from our point of view as 
theoreticians. The last review has been written in 1996 and much progress, on the theoretical, as well 
as, experimental side, has been made in the meantime. 

In |2U4j Austin gives a minireview of various indirect methods in nuclear astrophysics. He remarks 
that "it is a common perception that experimental nuclear astrophysics involves long measurements 
of small cross sections at lower and lower energies, so as to permit a reliable extrapolation to actual 
astrophysical energies. This perception is only partially correct. Recent developments, especially of 
radioactive beams, often permit one to obtain equivalent information with higher energy beams. The high 
energy experiments commonly yield higher event rates and sometimes yield information not available in 
the classical approach. " Coulomb breakup is a good example. 

In non-resonant charged particle reactions the energy dependence of the cross section is dominated 
by the penetration of the Coulomb barrier. This energy dependence is usually factored out by defining 
the astrophysical S-factor 

S{E CM ) = E CM (j(E CM ) exp(2vr?7), (113) 

where 77 is the Coulomb parameter Eq. (|13p. 

We want to give here some examples of Coulomb breakup experiments of astrophysical relevance. We 
also (briefly) explain their astrophysical relevance and the theoretical background of the experiments. 
In judging the conditions of an actual experiment, it is always useful to keep in mind the values of the 
basic parameters r], £, and x °f electromagnetic excitation, see Eqs. (fTU|) and (jUJ). In this way it is, 
e.g., quite clear that Coulomb dissociation of 8 B is a more favourable case than that of 16 0, with its 
comparatively high ^-values (see also Fig. 3 of |3*7j). 

7.1 Li Isotopes 
6 Li -> a + d 

The 6 Li Coulomb dissociation into a+d has been a test case of the Coulomb breakup method, see 
205 and fOJElCE]- These experiments were carried out with a 156 MeV 6 Li beam at the Karlsruhe 
Isochronous Cyclotron using the magnetic spectrograph "Little John" . This reaction is of astrophysical 
importance since the d(a, 7) 6 Li radiative capture is the only process by which 6 Li is produced in the 
standard primordial nucleosynthesis models. There has been renewed interest in 6 Li as a cosmological 
probe in recent years, mainly because the sensitivity of searches for 6 Li has been increased. 6 Li has 
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been found in metal-poor halo stars at a level exceeding even optimistic estimates of how much of it 
could have been made in the standard big bang nucleosynthesis. For more discussion on this, see 206J . 

A 6 Li-breakup experiment on 208 Pb at 60 MeV was performed at the Heidelberg tandem 207J. 
These authors found that the measured breakup cross section could not be related directly by first- 
order Coulomb excitation theory to the astrophysically relevant 4 He(a, d) 6 Li capture reaction. It would 
be interesting to redo the analysis with a modern reaction code and see the influence of higher order 
electromagnetic and nuclear effects. 

An experiment at the much higher beam energy of 150 AMeV than at Karlsruhe |205j or Heidelberg 
[207 is under way at GSI. Results from this experiments are eagerly awaited. 

For (N = Z)-nuclei electric dipole transitions obey the selection rule AT = 1. Since the two 
fragments d and a have isospin T = the electric dipole transition is forbidden by isospin selection 
rules. The contribution from the E2 transition should dominate. El transitions can still occur due to 
isospin violations. 

d(a, 7) 6 Li radiative capture experiments were also performed: Robertson et al. |20H] made measure- 
ments down to a relative energy of 1 MeV by detecting the recoil 6 Li ions. Cecil et al. |209j searched for 
the capture reaction at a relative energy of E cm = 53 keV; they obtained an upper limit on the reaction 
S'-factor of 2 x 10~ 7 MeVb. The direct capture into the 3 + resonance was measured also by Mohr et al. 

7 Li -> a + t 

The reaction t(a,7) 7 Li is also astrophysically relevant, and many experimental groups have measured 
this process in direct capture experiments. The nuclide 7 Li is produced in the early universe via 
this radiative capture reaction. The Coulomb dissociation method has also been applied under various 
conditions. We mention here the direct breakup of 70 MeV 7 Li scattered from a 120 Sn target |211j . It was 
found that the breakup is dominated by the Coulomb interaction for scattering angles inside the grazing 
angle. 7 Li breakup measurements on 197 Au at 54 MeV were performed in |212j . No straightforward 
analysis in terms of first order Coulomb breakup was found to be possible. 

A rather detailed measurement of the 7 Li — > a+t breakup was recently performed by Tokimoto et al. 
[73J with improved experimental techniques. This measurement was also accompanied by a theoretical 
analysis, where the time dependent Schrodinger equation was solved numerically, see [HI] and also the 
discussion in Sec. |3] In this way, valuable insight into the Coulomb dissociation process, the El and E2 
mixtures, its higher order effects (postacceleration) and the time-dependence of the tunneling process 
became possible. We refer the reader to this reference for further details. 

9 Li -> 8 Li + n 

The cross section for the radiative capture reaction 8 Li(n, r y) 9 Li was studied with the Coulomb disso- 
ciation of a 9 Li beam of 28.53 AMeV at MSU j2T5] . 

This reaction is of importance for the nucleosynthesis in inhomogeneous big bang models and in 
Type II supernovae. While the standard big bang nucleosynthesis ends with 7 Li, nucleosynthesis in 
neutron-rich regions could produce an observable amount of A > 8 nuclei. The 8 Li(n, 7) 9 Li reaction is 
in competition with the 8 Li(a, n) ll B reaction in determining the reaction path in the evolution network. 
This reaction may also play an important role in the explanation of the origin of light neutron-rich nuclei 
like 36 S, Ar, 46 Ca and 48 Ca. Recently, neutron-star mergers have been proposed as possible alternative 
sites for an r-process |214j . where similar reaction chains occur. 

This reaction is also of importance in order to determine the different primordial abundances of 
Li, Be, B, and C, in order to confine the possibility of inhomogeneous big bang nucleosynthesis. In 
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Figure 22: (a) The (n, 7) capture cross sections on 14 C are shown as a function of the cm. energy. 
The diamonds enclosed by two solid curves define the region of the experimental results. Plotted as 
1/3 of the actual values the dashed curve is a prediction of [220J, the dot-dashed curve is a prediction 
of |221j . (b) Blowup of the low energy region showing the same theoretical curves (not divided by 3), 
the experimental result of u 222 (solid point) and the lowest energy points (diamonds) of 182] with a 
kinematical fit (solid curve). This is Fig. 8 of 182 . 

such models the neutron-to-proton ratio is different in different regions in the early universe. In the 
neutron-rich region this reaction could therefore lead to a different bridging of the mass number A = 8. 

The half-life of 8 Li is less than a second (838ms), this probably makes a direct measurement of the 
capture cross section impossible. One has to rely on indirect methods, like the Coulomb dissociation of 
9 Li (with a half-life of 178 ms, but the Coulomb dissociation still works for such the 9 Li beam 

is produced in a fragmentation reaction). In this experiment |213j . only an upper limit could be put on 
the relevant cross section. Still, this is interesting since it directly rules out two theoretical calculations 
|215| I216j while the found upper limit is still consistent with the theoretical results of |217t [218J. A 
new attempt was recently made at MSU |219j to measure the 8 Li (7,n) 9 Li cross-section by the Coulomb 
dissociation of 9 Li with a more sensitive equipment than in |213j . It considerably lowers the upper 
bounds found there. We refer to this reference for details. 

n Li ^ 9 Li + In 

This reaction is mainly of interest for nuclear structure studies due to the Borromean nature of 11 Li. 
It is discussed in Sec. IB] above. Of general interest with respect to astrophysics is the question of the 
correlation of the two halo neutrons. This influences the rate of the two-neutron capture reaction. An 
understanding of this correlation is also of importance for other two-neutron capture reactions to be 
discussed in the Sec. 18.21 

7.2 14 C(n, 7 ) 15 C 

The nuclear structure aspect of the Coulomb dissociation of 15 C was touched upon in Sec. El where 
a measurement at GSI with a 606 AMeV 15 C beam is discussed. In this experiment, the B(E1)- 
distribution in the range of excitation energies from the threshold (E t hreshoid = 1.2184 MeV) up to 
about E = 8 MeV was studied. 

Now we are interested in thermonuclear energies, i.e., relative n+ 14 C energies from about 10 keV up 
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to about 300 keV (corresponding to temperatures between 0.1 — 3xl0 9 K). The Coulomb dissociation 
of 15 C was discussed in jTH] as a way to investigate the radiative capture reaction 14 C(n, 7) 15 C. This 
reaction is important in neutron-induced CNO cycles of stellar evolution phases beyond the main se- 
quence. It is also relevant in inhomogeneous big bang scenarios. A 15 C Coulomb dissociation experiment 
was performed at MSU |182j where details and further references can be found. The present status is 
not entirely conclusive, as one can see from Fig. E21 (Fig. 8 of [182 ). The excitation function for the 
14 C(n, 7) cross section was measured in |182j with an E/A = 35 MeV beam. This energy is much lower 
than the energy of 606 AMeV used at GSI. It is still high enough to break up the 15 C nucleus and 
excite the low energy continuum relevant for astrophysics. The deduced excitation function is shown 
as diamonds in Fig. l2*2T a). They are enclosed by the two solid curves, which define the region of the 
experimental results of the Coulomb dissociation experiment of |182j . Plotted as 1/3 of the actual 
values, the dashed curve is a theoretical prediction of Wiescher et al. [220 and the dot-dashed curve 
is a prediction of |221j . In Fig. l2*2T b[) a blowup of the low energy region is shown. The result from 
the direct capture experiment at 23 keV of the Karlsruhe group |222j (solid point) and the low energy 
points of the MSU experiment with a kinematical fit (solid curve) are shown. For further details see 

mi- 

The excitation function for the 14 C(n, 7)-reaction, as deduced from the Coulomb breakup, is con- 
sistent with an E}J 2 rise, as expected for p-wave capture. A fit gives <7 nj7 (23keV) = 3.2 ± 0.9/ib, about 
three times higher than the result of |222j . Also, the shape of the excitation curve has a peak at around 
E n = 200 keV. On the other hand, theoretical predictions by Wiescher et al. |220j and Descouvemont 
221 increase monotonically in this energy range. Another issue is the radiative capture contribution to 
the 5/2 + -excited state at 740 keV. This contribution is not included in the Coulomb breakup measure- 
ment. It is estimated to be small on theoretical grounds: the ^-dependence of the dipole transitions 
and the measured spectroscopic factors favour the ground state transition. We conclude that more work 
is necessary to resolve these questions. 

7.3 16 a J r 12 C Coulomb Dissociation 

The helium burning reaction 12 C(a, 7) 16 at thermonuclear energies is a key process for the evolution of 
massive stars and for the nucleosynthesis of 16 and heavier elements up to Fe. Despite its importance for 
nuclear astrophysics, the cross section for this reaction is still quite uncertain in the stellar energy domain 
(E cm ~ 300 keV). The reason is the extremely small value of this cross section and the superposition 
of El- and .E2-multipole contributions in the capture process. Several direct measurements of the 
radiative capture reaction, elastic a+ 12 C scattering and 16 N-decay studies have succeeded to determine 
reasonably well the .El-part of the astrophysical S-factor. However, the .E2-part, which is thought to 
be of about the same magnitude as the .El-part, is still quite uncertain. 

In this situation, Coulomb breakup seems to be an interesting method, since .E2-excitations are 
enhanced as compared to El, see Sec. especially Eq. (|2~3~j). However, several points render this method 
quite challenging. The a+ 12 C-breakup threshold E = 7.162 MeV is quite high, which means that the 
£- values, even at high 16 beam energies are not small (for the experiment mentioned below we get a 
value of about £ = 0.8). The relatively low flux of equivalent photons at these high energies leads to a 
serious competition with nuclear excitation effects. Also, there is El- and .E2-interference. Undaunted 
by all this, a breakup experiment was performed at GANIL at a projectile energy of 95 AMeV. Assuming 
only E2 Coulomb and AL = 2 nuclear excitation the measured breakup cross sections were compared 
to optical model calculations and an E2 5-factor was extracted in |223j . This analysis is based on 
theoretical calculations with the ECIS-code in 122 II and calculations of the fragment angular correlations 
in |225j . Further experimental studies were done recently at KVI , Groningen at 80 AMeV [226J. 
The large angular and momentum acceptances make the spectrometer BBS a very suitable device for 
such coincidence studies. A further experiment is planned at MSU which can use all the experience 
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gained from the previous experiments. The difficulties of the experiment were identified in the previous 
approaches and can hopefully be overcome. Much is at stake. As Wolfgang Pauli said (when he invented 
the neutrino): "Nur wer wagt, gewinnt!" ("Nothing ventured, nothing gained"). 

7.4 Some (j>, 7)- Reactions studied with Coulomb Dissociation 

The Coulomb breakup of 14 has been studied at RIKEN |227| and GANIL [I22B| in order to determine 
the S'-factor of the 13 N(p, 7) 14 capture reaction. It is astrophysically relevant for the hot CNO cycle. 
The 5-factor is dominated by the 1~ resonance in 14 which can be reached from the ground state in 
an El transition. Such a case is certainly simpler to study than a transition into the (flat) continuum. 
We refer to previous reviews [III EHj. We note that there is good agreement with the direct radiative 
capture measurement at Louvain-la Neuve [229 . This consistency can also be considered as a test 
of the Coulomb dissociation approach. It is also worth mentioning that as a by-product of the 14 
Coulomb dissociation approach the transition from the | ground state to the first excited \ + - state 
in 13 N was studied in |227j . The corresponding value of the radiative width T 7 was found to be in 
good agreement with values from a different previous experiment. Again, this can be considered as a 
confirmation of the Coulomb dissociation method. The accuracy of the RIKEN measurements [227 
was improved in new experiments at RIKEN |230| . There is agreement with the previous studies; the 
reliability of the Coulomb dissociation method has been tested to an accuracy of the order of 10 percent 
in the 12 C(p, 7) 13 N case. 

The 12 N — > n C + p breakup reaction was studied at GANIL |231| in order to study the n C(p, 7) 12 N 
radiative capture process. This reaction is relevant for the hot pp-chain. From the experimental breakup 
yield, the radiative width of the 1.19 MeV level in 12 N and the spectroscopic factor for the direct proton 
capture on 11 C have been extracted. The radiative width of the 1.19 MeV level is found to be smaller by 
more than one order of magnitude compared to a recent theoretical calculation but in rough agreement 
with an estimate by Wiescher et al. |232| . We refer to this reference for further details. The u C(p, 7) 12 N 
reaction was also studied at RIKEN by the Coulomb dissociation method |233| I234| . Their result is 
consistent with the GANIL measurement [231] , the accuracy of the value for the radiative width of the 
1.19 MeV level in 12 N is improved; for further details see these references. An experimental study of 
the 22 Mg(p, 7) 23 A1 reaction by the Coulomb dissociation method is also reported at RIKEN |235j . This 
reaction is relevant for the nucleosynthesis of 22 Na in Ne-rich novae. 

A very instructive test case would be 17 F-Coulomb dissociation. There are very accurate radiative 
capture experiments |236j and |237j to which one could compare possible Coulomb breakup measure- 
ments. The radiative capture reaction goes to the ground state and the | + excited state, with the 
Coulomb breakup method one can only study the ground state transition. There are data for the 
170 MeV 17 F dissociation into p+ 16 O 238J. They were analyzed theoretically by Esbensen and Bertsch 
[22J; substantial higher order effects were found at this energy. Higher 17 F beam energies would re- 
duce these effects and greatly facilitate the analysis. A recent work |239j found large nuclear effects 
in the dissociation of 17 F on 208 Pb at an energy of 65 AMeV. This would make an extraction of the 
electromagnetic transition matrixelement quite difficult. 

7.5 S B Coulomb Dissociation and the Solar Neutrino Problem 
Astrophysical motivation 

Two kinds of extraterrestrial neutrinos, solar neutrinos and supernovae neutrinos, have been observed 
and their measured fluxes, as well as, energy spectra, shed a new light on astrophysics and particle 
physics [240 . In this context, the 7 Be(p,7) 8 B radiative capture reaction is relevant for the solar neutrino 
problem. It determines the rate of production of 8 B which leads to the emission of high energy neutrinos. 
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The SNO experiment |241| has measured the v e flux in the charged current (CC) reaction 



d + v e ^p + p + e , 



(114) 



and also the neutral current (NC) reaction 



d + v. 



X 



p + n + v x 



(115) 



was recently measured (242[ 1248] . The CC reaction is only sensitive to electron neutrinos, whereas the 
NC reaction is equally sensitive to all neutrino flavours v x where x — e, fi, and r. 

The Kamiokande and SNO detectors also measure neutrino elastic scattering (ES) on the electrons: 



This reaction is sensitive to all kinds of neutrino flavours, mainly to electron neutrinos but also to some 
extent (to about 14%) to and v T via the neutral current interaction. Comparison of the two neutrino 
fluxes from SNO and Kamiokande and of the NC and CC results of SNO are a direct proof that there 
are neutrino oscillations. Some of the Kamiokande and SNO neutrinos come from flavours other than 
the electron flavour 3 . 

In this argument the astrophysical S-factor Syj does not enter, since one considers only the ratios of 
fluxes. However, this S'-factor is still of vital importance: it determines the absolute value of theoretical 
calculations of the solar 8 B v flux and therefore many efforts were undertaken to study the 7 Be(p, 7) 8 B 
reaction [213 12H]- 

In scenarios involving sterile neutrinos it is important to know the standard solar model 8 B neutrino 
flux in order to derive which part of the flux has oscillated into the sterile ones. We quote from 
247 : "Thus, the S B plays a crucial role in the interpretation of these experiments. Unfortunately, the 
predicted value of the 8 B flux normalization is quite uncertain, mainly due to poorly known nuclear cross 
sections at low energies. " 

For recent developments we refer to |248j . The situation is summarized in |249j : Improved ( S B) 
production rate predictions are very important for limiting the allowed neutrino mixing parameters, 
including possible contributions of sterile neutrinos. The astrophysical S factor 5*17(0) for this reaction 
must be known to ±5% in order that this uncertainty not be the dominant error in predictions of the 
solar v e flux. 

Direct measurements of the 7 Be(p, 7) 8 B-capture cross section 

There are direct reaction measurements, for the recent ones see Ref . |25U| I251| 1249} I252j . The target is 
the radioactive nucleus 7 Be and a major problem is the determination of the thickness of this target. 
The reaction has been measured at E c , m _ = 185.8,134.7 and 111.7 keV and the zero energy S'-factor 
inferred from these data is 18.5±2.4 eVb. This reaction has been studied recently in .253J with a 7 Be 
radioactive beam. In this way the target thickness problem is overcome, see also the discussion in 204J . 

There is an experiment by Junghans et al. from Seattle and TRIUMF [249] with a somewhat high 
value of the S'-factor. The cross section values of this paper are currently being revised, see Ref. 11 of 
[252 . A precision measurement at Ei a b = 991 keV with an implanted 7 Be target was recently reported 
in [232] . 

3 Oscillations of antineutrinos (coming from quite a few nuclear power reactors about 180 km from the detector) were 
recently also found by the KAMLAND experiment |244| . 



V x + e 



v x + e 



(116) 
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Coulomb dissociation of 8 B 

While there are certainly experimental problems with the direct capture experiments, there are (in 
addition to the specific experimental problems of the Coulomb dissociation method) problems of the 
theoretical analysis of the experiments in terms of the astrophysical S-factor. These theoretical problems 
can be solved in principle, as peripheral reactions are the ones best understood on a quantitative 
basis. They are mainly related to higher order electromagnetic effects (see also Sec. El above), nuclear 
effects (see Sec. EJ) and the mixture of different multipoles (especially El and E2). Yet, the Coulomb 
dissociation of 8 B offers some advantages: at the low binding energy of 8 B of 137 keV, the relevant 
photon fluxes are very high. Experiments under various conditions have been performed, which are 
affected in different ways by the different problems mentioned above. With all these problems identified 
a consistent picture emerges and an S'-factor is obtained in an entirely different way as compared to the 
direct measurements. This in itself is valuable. 

The story of the mid-nineties has been vividly told by Taube [254j and efforts to determine the 
astrophysical S-factor Sn by many different methods including Coulomb dissociation continue up to 
now. 

We now discuss the various experimental and theoretical contributions: Higher order perturbation 
theory using the semiclassical approximation was studied in ^3JE3- The semiclassical time-dependent 
Schrodinger equation was solved in [7UI73]. CDCC calculations were performed in [70], which become 
increasingly cumbersome for higher beam energies. The reactions at the different radioactive beam 
facilities have been done under different kinematical conditions, the highest 8 B beam energy is at GSI, 
254 AMeV [255J, at MSU it is about 40 to 83 AMeV |23il|25Z|. At RIKEN it is about 50 AMeV |25B| . 

The measurement at Notre Dame is at a very low 8 B beam energy of 26 MeV j259| . For such 
beam energies higher order effects are very important and it is not a simple task to extract a model 
independent S-factor. Nevertheless, these experiments are of great interest since they allow to test 
the dynamics of the Coulomb breakup process. With modern reaction codes (e.g. solving of the time- 
dependent Schrodinger equation [Bl^ 12*2*1 IT3 9] . see also Sec. EJ). One can well describe such reactions, 
however, the connection to the astrophysical S-factor is becoming somewhat indirect. 

One of the main issues is the contributions of the El and E2 multipoles and their separation. As 
can be shown from Eq. ()23|) the E2 photon number at intermediate and lower energies is enhanced, 
therefore they can contribute considerably to the Coulomb breakup cross section, whereas they are a 
small contribution to the photodissociation (or radiative capture) cross section. At RIKEN, one uses the 
different angular behaviour of the El and E2 multipolarities in order to disentangle the contributions, 
sec [260, 1261] I262| . The E2 contribution shows up predominantly at smaller impact parameters, i.e., at 
larger scattering angles. In Fig. |22]we show some theoretical calculations related to this point jT^j. The 
El contribution shows a (sin #) -2 -behaviour, whereas the E2 contribution is almost constant, down to a 
minimum angle & a diabatic- Nuclear diffraction effects are small. This is used in the RIKEN experiment to 
determine the E2 contribution. We show their result in Fig. These authors conclude that "although 
systematic effects could be large, the extracted E2 components appear to be quite small". 

In the MSU method, one uses the asymmetries, which are induced by the E1-E2 interference effects 
in order to sensitively determine the E2 S-factor [263. 12571 1264j . We show their longitudinal momentum 
distribution here in Fig. 

An intuitively appealing picture of the Coulomb dissociation process can be gained |257j by plotting 
b 2 Pdiss(b), where Pdi SS {b) is the Coulomb dissociation probability with various 7 Be scattering angle cuts. 
This is shown in Fig|2Hl one can clearly see that Coulomb dissociation essentially takes place outside 
the nuclear interaction region. 

In a particle-core model one can see qualitatively the relative importance of the various multipolar- 
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Figure 23: Coulomb dissociation cross section of 8 B scattered on 208 Pb as a function of the scattering 
angle for projectile energies of 46.5 AMeV (left) and 250 AMeV (right) and a 7 Be-p relative energy of 
0.3 MeV. First order results El (solid line), E2 (dashed line) and E1 + E2 excitation including nuclear 
diffraction (dotted line). (From Figs. 4 and 5 of Ref. |15j.) 




Figure 24: Observed cross sections eda/d8 are shown as a function of the outgoing 8 B* scattering angle 
for the three indicated relative energy bins. The y axis shows the product of the detection efficiency e 
and the differential cross section. The solid curve represents the best fit obtained with calculated El 
(dashed curve in (c)) and E2 (Coulomb and nuclear: dotted curve in (c)) amplitudes. For the first two 
energy bins shown in (a) and (b), the best fits result in pure El transitions. Dashed and dotted curves 
in (a) and (b) correspond to the results calculated with I = 1 and I = 2 components predicted by some 
theoretical models. This is Fig. 2 from |260j . where further details can be found. 
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Figure 25: The measured longitudinal momentum distributions of 7 Be fragments from the Coulomb 
dissociation of 44 AMeV 8 B on Pb with several maximum 7 Be scattering angle cuts are shown. The 
measurements are compared to a first order perturbation theory calculation convoluted with the ex- 
perimental resolution. The asymmetry due to the El-El interference is clearly seen. The overall 
normalization and the . )roduced from Fig. 8 of 

|257j . to which we refer Physical Society. 
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Figure 26: The product of the Coulomb dissociation probability P(b) and b 2 , where b is the impact 
parameter are shown as a function of b. Various 7 Be scattering angle cuts are applied, as indicated in 
the figures. Reproduced from Fig. 10 of |257j . Copyright (2001) by the American Physical Society. 
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Figure 27: Yields of breakup events (cross section times efficiency) plotted as a function of the relative 
energy. The solid and dashed histograms denote simulated E1+E2 and El yields respectively. Although 
the Ml contribution is reduced by about a factor of (f/c) 2 in the Coulomb dissociation the Ml peak 
(smeared out by the experimental resolution) is clearly seen in the GSI data. Reproduced from Fig. 1 
of |255j . Copyright (1999) by the American Physical Society. 



ities. The effective charges are given by 

Z$=*(_!^)\*pp-)\ (117) 
13 \m b + m c / \m b + m c / 

with the charge numbers Z b and Z c and the masses mj and m c of the fragments. We see that El charges 
are relatively small, E2 charges are large: for 8 B = 7 Be + p we have Z£ f = 3/8 andzg^ = 53/64. This 

can be compared to a neutron-core system: for 11 Be = 10 Be + n we find Z$ f = 4/11 and zfj f = 4/121. 
Furthermore, Ml transitions are suppressed by about a factor of (v/c) 2 . This Ml peak is nicely seen 
in the GSI experiment. It is somewhat smeared out due to the resolution in E re i in the p+ 7 Be invariant 
mass, see Fig. 

A new experiment on 8 B Coulombdissociation was recently completed at GSI [266 . They also 
pointed out that most extrapolations currently done are based on the cluster model of Descouvement 
267 , which does not seem to fully reproduce the energy dependence. This could partly explain the 
discrepancies of the different experiments done at different relative energies. This idea is pursued further 
in |268j taking into account a full dynamic reaction calculation, as well as, an improved potential model 
for 8 B. A recent minireview can also be found in j269j . where also a future CERN/ISOLDE measurement 
using a radioactive 7 Be beam is briefly discussed. 

The experimental results was summarized in Fig 19 of |257j . We have updated this figure with 
the latest results and show it in Fig. |2H Altogether it is quite remarkable that completely different 
experimental methods with possibly different systematic errors lead to results that are quite consistent. 



8 Possible New Applications of Coulomb Dissociation for Nu- 
clear Astrophysics 

Radioactive beam facilities all over the world have considerably widened the scope of nuclear physics 
in the past years. In the previous sections we described the applications of Coulomb excitation and 
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Figure 28: The zero-energy astrophysical S'-factor S 17 for the 7 Be(p, 7) 8 B reaction is shown from selected 
direct and indirect methods. The open circle denote results obtained from direct measurements, the full 
circle those of Coulomb breakup experiments. The star is the value given in |265| as proposed reference 
value. Adapted from Fig. 19 of |257j . 

dissociation to nuclear structure and astrophysics. For example halo nuclei and new regions of defor- 
mation could be studied in detail for light exotic nuclei. With the new forthcoming generation of RIA's 
also medium and heavy nuclei will be produced with sufficient intensities, with further opportunities. 
We mention the project in Europe where a substantial part of the program is in Rare Isotope 
Beams research. In the US the scientific opportunities with fast fragmentation beams from the Rare 
Isotope Accelerator are intensively investigated see jlj, and also |270j . In Japan, at RIKEN, Wako the 
construction of a Radio Isotope Beam Factory (RIBF) is underway, see e.g. [0]. 

A status report is also provided in |271j . The astrophysical scenario is presented by Kappeler et al. 
in |272j . In the following we indicate some of the future opportunities in this field related to nuclear 
astrophysics, especially to the r- and rp-processes. 

8.1 Electromagnetic Properties of r -Process Nuclei 

Nucleosynthesis beyond the iron peak proceeds mainly by the r- and s-processes (rapid and slow neutron 
capture). This is widely discussed in the literature, see e.g. |273t 1274] or also [275 . Concise executive 
summaries are provided in |3| (p. 96f) or j3] (p. 54). 

In order to assess the r-process path, it is important to know the nuclear properties like /3-decay 
half-lives and neutron binding energies. In the waiting point approximation |273| I274j an (n,7)- and 
(7,n)-equilibrium is assumed in an isotopic chain. As the nuclei inside the isotopic chain are assumed 
to be in a thermal equilibrium, only binding energies and /5-decay halflives are needed. 

In general the waiting point approximation should be replaced by a dynamic r-process flow cal- 
culation, taking into account (n/y), (7,n) cross sections and /3-decay rates, as well as, time- varying 
temperature and neutron density for the astrophysical scenario. In this case the knowledge of (n,7) 
cross sections is of importance. 

To establish the quantitative details of the r-process, accurate energy-averaged neutron-capture 
cross sections are needed. Such data provide information on the mechanism of the neutron-capture 
process and time scales, as well as, temperatures involved in the process. The data should also shed 
light on neutron sources, required neutron fluxes and possible sites of the processes (see Ref. |273p . In 
this situation, it is very instructive to look at the model studies of [276 and the scenario described in 
|277j . The dependence of direct neutron capture on nuclear structure models was investigated in [276J . 
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Figure 29: Dependence of level energies on mass number for even-odd Sn isotopes calculated in the 
Hartree-Fock-Bogolyubov model (HFB) (upper panel), in the Relativistic Mean Field approach (RMF) 
(middle panel) and in the Finite Range Droplet Model (FRDM) (lower panel). The l/2~state (open 
circles), the 3/2~ state (triangles) and the calculated neutron separation energy (full circles) are shown. 
The lines are drawn to guide the eye. These are adapted from the Figs. 10-12 of [276J. 



57 




Figure 30: The measured and calculated abundances of nuclei are compared with each other. The dots 
show the measured abundance which are compared to the one of a calculation including a low energy 
dipole mode and one with a normal GDR mode. The presence of a low energy GDR leads to higher 
abundances for larger A (lighter curve), whereas the one with a standard GDR mode has a higher 
abundance at lower A (darker curve) This is Fig. 24 of 0], which is an adaption from |277j . 

The investigated models yield capture cross sections sometimes differing by orders of magnitude. The 
A-dependence of level energies and neutron separation energies is illustrated by comparing Figs 10-12 
of |2ZE!, shown in Fig. 

The dependence of the level energies and neutron separation energies on various theoretical models 
(HFB, RMFT and FRDM) can strikingly be noticed there. This may also lead to differences in the 
predicted astrophysical r-process paths. Because of the low level densities, the compound nucleus model 
may not be applicable. 

The dependence of the r-process abundance on nuclear properties and astrophysical conditions is 
investigated in |277j . The influence of the normal GDR and a GDR, which includes low energy dipole 
contributions, is strikingly summarized in Fig. EH Specifically the possible existence of a low-enery El 
pigmy resonance is studied. There can also be a problem with the application of the statistical (Hauser- 
Feshbach) model to neutron rich nuclei: since the number of available resonance states may not be large 
enough for the application of this model, there could be a possible overestimate of the statistical (Hauser- 
Feshbach) predictions for resonance-deficient nuclei. The influence of nuclear structure properties on 
the r-process abundance distributions is shown in Fig. 6 of |277j for the astrophysical scenarios T = 
1.0 x 10 9 K, N n = 10 20 cm" 3 , r = 2.4 s and T = 1.5 x 10 9 K, N n = 10 28 cm" 3 and r = 0.3 s, where N n 
is the neutron density and r the irradiation time. We recall that the astrophysical site of the r-process 
is still under debate. From this study it can be concluded that |277j "there is an urgent need to carry 
on investigating theoretically as well as experimentally, the radiative neutron captures by exotic nuclei 
in order to improve our understanding of the r-process nucleosynthesis. " 

In such a situation, the Coulomb dissociation can be a very useful tool to obtain information on 
(n,7)-reaction cross sections on unstable nuclei, where direct measurements cannot be done. Of course, 
one cannot and need not study the capture cross section on all the nuclei involved; there will be some 
key reactions of nuclei close to magic numbers. It was proposed in |278j to use the Coulomb dissociation 
method to obtain information about (n,7) reaction cross sections, using nuclei like 124 Mo, 126 Ru, 128 Pd 
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and 130 Cd as projectiles. 

Since the flux of equivalent photons has essentially al/w dependence, low neutron thresholds are 
favourable for the Coulomb dissociation method. The optimum choice of beam energy is an important 
issue. It will depend essentially on the actual neutron binding energy. There is a factor of {c/v) 2 on 
the photon flux (see Eq. (|20p) favoring the low beam velocities. On the other hand, the beam energy 
should not be so low that higher order effects become important (the strength paramer x has an In- 
dependence, see Eq. ((12} ) or that the adiabaticity criterion is not fulfilled (see Eq. (|T0|) ). Note that 
only information about the (n.,7) capture reaction to the ground state is possible with the Coulomb 
dissociation method. The situation is reminiscent of the loosely bound neutron-rich light nuclei, like 
11 Be, n Li and 19 C, which were discussed in Sec. We recall that it was possible to determine the 
neutron binding energy of 19 C from the shape of the low-lying El strength distribution, see (HHl or 
Sec. El The relevant level densities are low, which is also favourable for the Coulomb dissociation 
method. To first order Coulomb dissociation probabilities are independent of the mass of the core (the 
recoil is proportional to 1/A, but the Coulomb force is proportional to Z, see Eqs. ()I2J) and (jlHI)), so 
one can expect that the Coulomb dissociation method works essentially in a quite similar way as for 
the light ions like u Be, or 15 . 17 > 19 C. 

In [T3| the I st and 2 nd order Coulomb excitation amplitudes are given analytically in a zero range 
model for the neutron-core interaction (see also Sec. |HJ). We propose to use this handy formalism to 
assess, how far one can go down in beam energy and still obtain meaningful results with the Coulomb 
dissociation method, i.e., where the 1 st order amplitude can still be extracted experimentally without 
being too much disturbed by corrections due to higher orders. For radioactive beam facilities, like ISOL 
or SPIRAL, the maximum beam energy is an important question for possible Coulomb dissociation 
experiments. (For Coulomb dissociation with two charged fragments in the final state, like in the 8 B 
— > 7 Be+p experiment with a 26 MeV 8 B beam |259j it seems to be impossible to obtain such a simple 
analytical formula and one should resort to the more involved approaches mentioned in Sec. El) 

There is also another aspect: the radiative neutron capture by neutron rich nuclei at astrophysical 
energies is given by the low energy tail of the El strength. Standard values for the Lorentz parameters 
of the GDR are normally used. The possible existence of a low energy El pigmy resonance can strongly 
influence this picture, see |277j . Coulomb dissociation is the tool to study this question experimentally; 
we refer to our discussion of the light ions in Sec. El Especially, at the GSI the dipole response of the 
oxygen isotopes has been studied |189j . It will be very interesting to see the corresponding results for 
the (medium) heavy nuclei. 

With the new radioactive beam facilities (either fragment separator or ISOL-type facilities) some 
of the nuclei far off the valley of stability, which are relevant for the r-process, can be produced. An 
impression of the future possibilities can be obtained, e.g., from Fig. 2 p. 10 of j3]. In the chart of 
nuclides the production rates predicted for RIA can be seen, along with a line which indicates the 
r-process path. The rates for the future accelerator at the GSI are given in Fig 1.17, p. 106 and on 
p. 156 in 

The estimated minimum intensities of high-energy secondary high energy beams necessary for a 
certain type of reaction studies are given in Table 1.2, p. 125 of [3]. For Coulomb breakup a rate of the 
order of 10 3 ions/s is given as a rough general estimate. 

8.2 Two-Particle Capture Reactions and Cross Sections Relevant for the rp-Process 

A new field of application of the Coulomb dissociation method can be two particle capture reactions of 
the type 

A + P! +p 2 -> 5 + 7. (118) 

The most famous reaction of this type is certainly the triple a capture reaction leading, via an 8 Be 
unstable state, to a resonance in 12 C, which was predicted by Hoyle |279j . Evidently, such a type of 
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Figure 31: The rp-process path, including 2p-captures, for temperatures of 1.9xl0 9 K and densities of 
10 6 gcm~ 3 |285j . Also shown are stable nuclei and the position of the proton-drip line. This is Fig. 1.14 
of0. 

reaction cannot be studied in a direct way in the laboratory. Sometimes this is not necessary, when the 
relevant information about resonances involved can be obtained by other means (transfer reactions, etc.), 
like in the triple a-process. Another way to investigate such processes can be the Coulomb dissociation 
method, where the time reversed process 7 + -B(gs) — > A + pi + p 2 is studied using equivalent photons. 
Let us give some examples: 

Two- neutron capture reactions in supernovae neutrino bubbles are studied in Ref. |280| . In the case 
of a high neutron abundance, a sequence of two-neutron capture reactions, 4 He(2n,7) 6 He(2n,7) 8 He can 
bridge the A = 5 and 8 gaps. The 6 He and 8 He nuclei may be formed preferentially by two-step resonant 
processes through their broad 2 + first excited states j2£0j. Dedicated Coulomb dissociation experiments 
can be useful. The 6 He Coulomb dissociation is measured in |18Sj. Astrophysical aspects related to the 
two-neutron capture on 4 He are also discussed there. One can envisage two mechanisms: a nonresonant 
two-neutron capture and a capture with an unstable 5 He resonant intermediate state |281j . A dedicated 
Coulomb dissociation experiment could shed light on the question which mechanism is dominant. 

Another key reaction can be the 4 He(an,7) 9 Be reaction |28(Jj . The 9 Be(7,n) reaction has been 
studied directly (see Ref. |282j ) and the low energy si resonance is clearly established. Despite this, a 
9 Be Coulomb dissociation experiment could be rewarding (see also Ref. |178p . Other useful information 
is obtained from (e, e') and (p, p') reactions on 9 Be [283.. Recently the photodissociation of 9 Be was 
studied by |284| 203 , as was already mentioned. 

Another important application is for reactions relevant for the rp-process: neutron-deficient nuclei 
close to the proton drip line play an important role in a variety of astrophysical scenarios such as Nova 
explosions, X-ray bursts, or X-ray pulsars. In those scenarios, hydrogen is burnt via a sequence of 
rapid proton captures and /3 + -decays close to the proton drip line (called the rp-process). The impact 
and perspectives of radioactive beam experiments for the rp-process was recently discussed by Wiescher 
and Schatz |286j . It is noted there that (p, 7) capture experiments become increasingly difficult with 
increasing charge Z, and Coulomb breakup, Coulomb excitation and particle transfer reactions in inverse 
kinematics become promising indirect methods to access the relevant information. We recall that (p, 7) 
reactions on light nuclei were studied by the Coulomb dissociation of 8 B, 12 N, and 14 radioactive 
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beams, see Sec. |7J For further discussion we refer to [3 j p. 98ff and especially Fig. 1.14, shown here as 
Fig.EU 

In the rp-process, two-proton capture reactions can bridge the waiting points |287| 12881 l^85j . From 
the 15 0(2p,7) 17 Ne, 18 Ne(2p, 7) 20 Mg and 38 Ca(2p, 7) 40 Ti reactions considered in Ref. [2HH], the latter 
can act as an efficient reaction link at conditions typical for X-ray bursts on neutron stars. A 40 Ti — > 
p+p+ 38 Ca Coulomb dissociation experiment should be feasible. The decay with two protons is expected 
to be sequential rather than correlated (" 2 He" -emission, 2p-radioactivity). The relevant resonances are 
listed in Table XII of Ref. |288j . In Ref. |285j it is found that in X-ray bursts 2p-capture reactions 
accelerate the reaction flow into the Z > 36 region considerably. In Table 1 of Ref. [285] nuclei, on 
which 2p-capture reactions may occur, are listed; the final nuclei are 68 Se, 72 Kr, 76 Sr, 80 Zr, 84 Mo, 88 Ru, 
92 Pd and 96 Cd (see also Fig. 8 of Ref. j2HZj)- It is proposed to study the Coulomb dissociation of these 
nuclei in order to obtain more direct insight into the 2p-capture process. 



9 Conclusions 

Peripheral collisions of medium and high energy nuclei (stable or radioactive) passing each other at 
distances beyond nuclear contact and thus dominated by electromagnetic interactions are important 
tools of nuclear physics research. The intense source of quasi-real (or equivalent) photons present 
in such collisions has opened a wide horizon of related problems and new experimental possibilities 
especially for the present and forthcoming radioactive beam facilities to investigate efficiently photo- 
interactions with nuclei (single- and multiphoton excitations and electromagnetic dissociation). We 
have described the basic points of the theory, which rests on very solid grounds: quantum mechanics 
and the electromagnetic interaction. Problems in the technical description of the process, problems 
due to nuclear interactions between the ions are identified and can be considered as well understood. 
Modern computational methods are of great help. 

It has sometimes been said about the Coulomb dissociation method, see e.g. [249 , that "it is difficult 
to determine all of their important systematic errors". Certainly, there are systematic effects which one 
has to take carefully into account. This is the task of nuclear reaction theory which is described in this 
review. We hope that it has become clear that these effects are well identified and can be calculated. 
After all, it is essentially the electromagnetic interaction which enters, and there is QED, the best theory 
we have. We hope that we have shown here that the view expressed in [249J is overly pessimistic. 

Certainly, there are also experimental problems: one of them is the accurate determination of the 
relative energy of the fragments. Due to the influence of the Coulomb barrier, the astrophysical ^-factors 
depend very sensitively on this quantity. However, we have to leave such questions to our experimental 
colleagues. 

After these theoretical considerations we discussed experimental results in the field of nuclear struc- 
ture and nuclear astrophysics. Due to the good theoretical understanding of the electromagnetic exci- 
tation, unambiguous conclusions have been drawn from them. Let us mention here the discovery of low 
lying El strength in neutron- rich nuclei and the determination of astrophysical S-factors of radiative 
capture processes, like 7 Be(p,7) 8 B. 

Fast moving heavy nuclei are a bright source of equivalent photons, and we expect a bright future 
of this subject, especially at the future rare isotope facilities. Coulomb excitation and dissociation is 
an opportunity to study the interaction of exotic nuclei with photons. It is unique, until electron-ion 
colliders will become available as a complementary tool. 
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A The Electromagnetic Interaction: the Condition of No 
Nuclear Contact and Expansion into Electromagnetic Mul- 
tipoles 

Although the following is quite familiar, it seems of interest to deal with some basic physical points 
which arise in the special case of electromagnetic excitation in nucleus-nucleus collisions. We especially 
wish to clarify the difference of it to the treatment of electron scattering, which is perhaps even more 
familiar to most readers. 

As was already emphasized, in hadron-hadron collisions the important condition for the dominance 
of the electromagnetic interaction is r < r p , i.e., there are two spatially separated charge distributions, 
see Fig. [IJ In this way strong interactions are avoided, there is only the electromagnetic interaction 
between the ions and this simple reaction mechanism is not spoiled by nuclear effects. It is useful 
to take the photon propagator in coordinate space f and only Fourier transformed with respect to t, 
i.e. we consider D^if, u); in this representation the condition of no nuclear contact (r < r p ) can be 
implemented easily. 

The condition r < r p is used in the semiclassical approach here, where the projectile on the orbit 
r p (t) does not overlap with the target. In a distorted wave Born approximation (DWBA) with, e.g., 
only Coulomb waves (CWBA) the no nuclear contact condition r < r p is also useful: for low collision 
energies, where the classical turning point is outside the target nucleus, the Coulomb wave function is 
very small inside the target nucleus and its contribution can be neglected. 

As we will see below, the information about the hadronic system to be studied can be expressed in 
terms of electromagnetic matrix-elements at the photon point. 

The interaction of two currents is a standard problem in electrodynamics and QED. The electro- 
magnetic field is described by the vector potential and its interaction as M int = J d 4 x j ■ A. We 
define the scalar product of two four-vectors as a ■ b = g^a^by = = a b — a ■ b, where we use 
the metric g = diagonal(l, —1, —1, —1) throughout this article. In our case the vector potential A^ is 
generated by the moving charge distribution J M . It is given by 




(119) 



The Fourier transform of the current distribution is defined as 




(120) 




with the notation d A k = dkod?k and k 



(ko = u,k). We have current conservation, which is expressed 



as 



Q^f = djo/dxo - V • j = 0, 



(122) 
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or k ■ j = in momentum space. 

The electromagnetic interaction between two current distributions j\ and j 2 can be written as 



M int = J J d 4 xd A x'j llx (x)D^(x - x')j 2v {x'). (123) 

Alternatively it is given in terms of the Fourier transform of the propagator 

D^k) = J d A xexp(-ikx)D uu (x) (124) 

(125) 

and the currents (see above) as 

M int = j^Jd A k] lfl (-k)D uu (k)] 2u (k). (126) 

In the Born approximation which is often used, e.g., in lepton scattering, see Fig. EJ Eq. (|126|) 
simplifies as we can use plane waves for initial and final states. We consider the reaction 1 + 2 — ■> V + 2', 
with the four- momenta pi+p 2 = v'\ +v'i- The momentum transfer q is given by q = p\ — p[ = —p 2 +p 2 . 
The transition currents can be written as ji u (x) = a M exp(i(pi — p[)x) and j 2fJi = b u exp(i(p 2 — p' 2 )x') 
where and 6 M are independent of x, x'. With this we can write Mi n t as 

M^t = (2vr) 4 5( Pl + p 2 - p[ - p' 2 )a n D nu (q)b u . (127) 

The fundamental property of the electromagnetic interaction is gauge invariance: We can add to 
any vector potential a gradient term of the form d^xix)- Due to current conservation d^j^ = this 
leaves the electromagnetic interaction M^t invariant. Various gauges exist also for D nu , corresponding 
to various gauges of A u . E.g., [501 gi ves a ver Y instructive discussion on the application of different 
gauges, see also |289j p. 109ff. For interesting historical remarks (e.g., that the Lorentz condition is 
originally due to Lorenz) see [290 . A rather general covariant form of the propagator is given by the 
linear lambda gauge 

, x 47r / ^k n k v \ , . 

D ^ k) = WTTe{ 9 ^ + x ^)> (128) 

where A is an arbitrary parameter. 

In the Coulomb gauge (characterized by divA = 0) one has explicitly jUHI 

47T 



00 



D 

k 2 

D 0l = Ao = 



Att / kiki\ 

D " = ~w[ s -'-w)- < 129) 

(with k 2 = u) 2 — k 2 and u = k ) where i, I — 1,2,3 denote the spatial indices. (We leave out the +ie for 
simplicity.) In this gauge we have 

a M k)K = ^\=^- S -^ + ^^). (130) 



Using current conservation we can rewrite this as 

= 47T 



a ■ b a b ( _ 



k 2 k 2 \ k 2 
k 2 V^o-a-6>! 



An / -\ An 

[a b -a-bj = —a ■ b. (131) 
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Due to gauge invariance one obtains the same result in the Coulomb gauge, as well as in the Feynman 
gauge (A = 0), Lorentz gauge (A = — 1) or in any other gauge, as it should be. 

In the general case, which is relevant for electromagnetic excitation with nuclei, the transition 
currents are not associated with a definite momentum transfer q. There is an integration over the 
momentum transfer k, see Eq. ()126j) above. In the static limit (no retardation, corresponding to 
k = uj = 0) one gets 



M, 



int 




ii,o(-fc)j2,o(fc) -ji(-fc) • 32{k ) 

k 2 

Pi{r 1 )p 2 {r 2 ) -Ji(ri) ; J2(r 2 )/c 2 
In - r 2 | 



(132) 



which corresponds to Eq. (0) in the main text, the starting point for the nonrelativistic Coulomb 
excitation in [T], see also p. 84 of [2*T] . 

Now we want to indicate how to express the interaction in terms of the electromagnetic multipole 
matrix-elements of the charge- current distribution. (The details can be found in the textbooks, see, 
e.g., jZIJ-) This is possible for non-overlapping charge distributions. We only deal with the monopole- 
multipole part of the interaction, see Fig. ^ Then we can assume that the charge Z p is a point charge. 
The static case is well known, the relativistic case can be handled as well. 

In [2T] the Lorentz condition is used for the propagator. See especially Sec. 6.1 Eq. (8) and (22) and 
Sec. 6.2 in that reference. We do not give here all the details of the multipole expansion for interacting 
charges and currents, see, e.g., [2TJ[3UH291J. Let us emphasize the main points related to the question of 
penetrating (like in the plane wave case) versus spatially separated charge distributions. The propagator 
("Green's function") is expanded into multipoles in analogy to Eq. (plj (see also, e.g., f2T|. Eq. (23)) 4 

expMf-r;|) = 4 ^^ jLK<)/l (i) K>)yZM( ^ <)riM( ^ >)) (133) 

l r r P\ LM 

where r< and r > are the smaller and larger of r p and r respectively. In the case of spatially separated 
charge distributions we always have r< = r and r > = r p . 

This is a generalization of the multipole expansion of the static Coulomb interaction 5 , 

T^FT = ^-7-^lYL*M(r<)YLM(r>), (134) 

which can be obtained from Eq. ()133j) by performing the limit uo — > and using the expansion 

^)=(^H)h> and ^) = -*^z# ( 135 ) 



4 We follow here the convention of 92 in the definition of the spherical Bessel function, which is {z) — jL{z)+iyL(z) 1 
where jh and jjl arc the (regular) spherical Bessel and the (irregular) Neumann function. This is different than the one 
used in [23 • 

5 For L = we have a term proportional to l/r>. This term gives a contribution to nuclear excitation only if r > r p . 
Monopole transitions can occur, e.g., in the case of internal conversion, due to the penetration of the electron charge 
cloud into the nucleus, see, e.g., |292| . In this case we have r > r e , where r e and r denote the electronic and nuclear 
variables, respectively. Usually, the electron cloud is outside of the nucleus, i.e. r < r e and the probability of internal 
conversion is proportional to the B(tt, A)-value of the corresponding 7-transition. This is in close analogy to the Coulomb 
excitation discussed here. 
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for small values of the argument x = ujt. 

On the other hand one can Fourier transform the (retarded) interaction into momentum space by 
means of the Bethe integral 

exp{iuj\r-f p \) _ 1 f ^ expjik ■ (r- f p )) 



\r-r p \ 2n 2 J k 2 - cu 2 - ie 

with e being an infinitesimally small positive number. Using the plane-wave expansion of exp(ik-f) into 
multipoles, integrating over the angular part of k, and making use of the orthogonality of the spherical 
harmonics we obtain the formula 



CO 



<f )] ; W = ^<<"<)k£W (137) 

K UJ %a A 



In the static limit uj = one obtains 

7r r 



OO rrr r^L 







dkj L (kr)j L (kr p ) = (138) 



In [2] the Coulomb gauge is used to derive the full expression for first order Coulomb excitation, see 
their Sec. II B.l. The result in this approach is the same as in [21], where the Lorentz gauge is used; this 
must be so due to gauge invariance. It is important to note that the electromagnetic matrix-elements 
at the photon point (see Eqs. (|141|) and (|142j) below) enter in the expression for the Coulomb excitation 
amplitude. In the approach using the Coulomb gauge one also needs the equation (Eq. (2.B.11) in [2]) 

dh k 2 -u 2 -ie = to 3L( > Ur < )hL {Ur>) ~ 2(2L + 1V^ ( 139 ) 
which can either be found in |293j or follows from the identity 

k 2 1 1 



o 



uj 2 (k 2 — uj 2 ) uj 2 k 2 — uj 2 ' 



(140) 



together with the two expressions Eqs. ()137|) and Eqs. ()138|) above. 

These expressions are now used in the calculation of the matrixelements for electromagnetic exci- 
tation, either in a semiclassical or in a quantum-mechanical framework (for all the necessary details 
on the multipole expansion of the electromagnetic field and on vector spherical harmonics see, e.g., 
[2*1 |2*H H7T| ). A good discussion of vector spherical harmonics and electromagnetic interactions is also 
given in |294j . One clearly sees how the separation into the r and r p parts comes about when we can use 
the no penetration condition r< = r and r> = r p . The expressions factorize naturally into an orbital 
part (associated with (ur p )) and an electromagnetic matrixelement (associated with the ji(ur)). 

The Coulomb excitation amplitude is expressed in terms of the electromagnetic multipole moments, 
which are defined as (see, e.g., [2]) 

M(EXfx,q) = J+lVfi J d 3 rj- V x L [j^Y^f)) (141) 

M{MXfjL,q) = -i ^+y J d 3 rj(7) ■ L [ji(qr)Y lm (r)) , (142) 

with L = -irxV (see Eqs. (II.1.3), (II.1.4) of P). 

One can see that the integration over the momenta k (see Eqs. ()137j) and ([138") ) conspires in such a 
way that only the electromagnetic matrixelements at the photon point k = uj (or the long-wave-length 
expression, respectively) appear. 
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An even more intriguing version of an addition theorem to expand the electromagnetic field into 
multipoles and into r< and r > parts is used by Winther and Alder [21] to handle the semiclassic 
relativistic straight line case. In this case the electromagnetic potential is described by the Lienard- 
Wiechert potential 

<f>(r',t)= . Zp€1 A(r',t) = -<f>(?,t), (143) 

^(b-x') 2 + y' 2 + 1 2 {z' -vt) 2 c 

which can be Fourier transformed with respect to t as 

2Z n e ( ,uj ,\ l uj 



Mjj oj) = —t- exp i-z' K — Mb - x') 2 + y' 2 ). (144) 
v V v J \wy v / 



One expands this expression into multipole components 



r 



\M = T,Wx,(r',uj)Y^(?). (145) 



By the use of the Graf addition theorem (see Eq. (9.1.79) of |52|) an analytic expression for 

is obtained. Again a complete separation of the excitation probability in electromagnetic multipole 

moments and quantities which describe the projectile motion is reached, see Eq. (123)1 of Sec. [2j 

This is to be contrasted to the plane wave case (e.g. commonly used in inelastic electron scattering). 
In this case there is a definite (space-like) momentum transfer q = p% — p[, see Fig. El By varying q one 
can probe the structure of an object in a way not possible with real photons. An important example is 
deep inelastic electron scattering. 

For small values of — q 2 we can think of the exchanged photon as quasireal. For small energy loss and 
small angle scattering we have the kinematical relation — q 2 = q\ + {ujji^v )) 2 . As the main contribution 
to the total cross section comes from small values of — q 2 (where the photon propagator becomes very 
large), this cross section is dominated by values of — q 2 ~ (uj/^v)) 2 , which, especially for large values 
of 7, can be quite small and therefore q 2 ~ (quasireal) photons dominate in this case as well. 

This kind of equivalent photon approximation is different from the one discussed mainly in this 
review: here we have used the semiclassical approximation, as it is appropriate for the nucleus-nucleus 
collisions. This leads to impact parameter dependent equivalent photon spectra. These spectra can be 
integrated over the impact parameters b > R m in- Still, the plane- wave and the semiclassical approach 
have some qualitative features in common, like the logarithmic rise of the photon number with the 
Lorentz factor 7. As the name implies, the EPA is an approximation, with certain ranges of validity, 
different for the semiclassical and plane wave approaches, see, e.g., |295j and |296j . In |297j experimental 
data for electromagnetic excitation with heavy ions were analysed using the EPA spectrum of [296J. 
In contrast to semiclassical calculations, systematically lower cross sections are obtained that cannot 
reproduce the experimental results. 

How the EPA works and where its limits are, is explicitly studied, e.g., in the case of H° production, 
298 , 299| l3UUj : we can only refer to these papers here. The PWBA in the Coulomb gauge is also studied 
very lucidly in a classic paper by Fano |3(Jlj : a more recent work is |302 . 
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